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ABSTRACT 


Subramaniam, Lavanya, <Bachelors of Engineering, College of Engineering, Anna 
University, India, May 1999>, Partitioning a Non-simple Polygon into Simple Polygons, 
Chair of Committee: Dr. Thomas F. Hain. 


While there are several algorithms to partition simple polygons into triangles, trapezoids, 
or monotone polygons, they do not handle self-intersecting (complex) polygons. In this 
work, we describe the theory, implementation, and testing of a new object-precision 
algorithm to partition a complex polygon into a set of simple polygons (possibly having 
holes) using the non-zero winding number rule to determine the “inside”’ness of points in 
the polygon. The outline for the algorithm implemented and tested in this research was 
developed by Dr. Thomas Hain. This algorithm can be used as a preprocessing step to 
convert complex polygons into a set of simple polygons, which are then amenable to 
being processed by the algorithms mentioned above. We show that this algorithm yields 
better high resolution performance than the image-precision scan line algorithm[6]—the 


method currently used in practice. 


vill 


Chapter 1 
INTRODUCTION 


1.1 Motivation for polygon partitioning 


Efficient polygon partitioning has been one of the most outstanding open 
problems in 2-dimensional computational geometry, and is usually confined to the 
domain of simple polygons. Simple polygons are ones that have no self-intersecting 
edges, but that can nevertheless be quite convoluted. To solve a variety of problems, the 
target polygon is typically partitioned into a disjoint set of covering polygons of restricted 
form, such as triangles, trapezoids, convex polygons, monotonic polygons, or polygons 
having other desirable characteristics. 

Thus, polygon partitioning is one of the pre-processing steps for many 
applications in computational geometry and hence in computer graphics. Some of the 
applications requiring such polygon partitioning are: 

e Finding the area of a simple non-convex polygon. 

e Discovering whether a point is inside a simple polygon. 
e Finding the center of gravity of a simple polygon. 

e Filling polygons. 


e Discovering whether a polygon is simple (has no self-intersecting edges). 


Some real-time applications have triangulation as one of their major components. 
For example, the Art Gallery problem, which deals with monitoring a convex region, 
involves triangulation to minimize the number of cameras required to cover the whole 


guarded area. 


Although the partitioning problem seems simple, the solutions tend to be quite 
complicated. These problems are often complicated and there has been a lot of research 
[2], [9], [17] that has focused on achieving high efficiency, or low running-times. 


1.2 The motivation for the current work 


The goal of this research is to offer an efficient algorithm to partition complex 
polygon into simple polygons so that the domain of such sophisticated algorithms already 
developed, can be extended to include complex polygons as shown in Figure 1. In 
addition, there are polygon-processing algorithms, which result in complex polygons, 


which can be further processed after the application of this algorithm. 
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Figure 1 Domain of current work 
In a larger context, we have taken as our domain of interest the printer industry. 


Here, high-level page objects (such as complex polygons) to be rendered are described in 
a page description language (PDL), such as PostScript®, HP PCL, or HP-GL. These 
objects are given to a graphics library which converts them into low-level graphical 
primitives (such as trapezoids) that are then stored in a display list. The objects in the 
display list can be quickly scan-converted (turned into pixels) into a printed image. The 
display list represents a semi-digested and compressed version of the page description, 
and can be efficiently scan-converted multiple times for multi-page, collated printing. 
The algorithm described here would reside in the graphics library, together with an 
already existing trapezoidation algorithm to convert simple polygons into lists of 


trapezoids. 
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Figure 2 Printer graphics [3] 


1.3. Definitions 


A simple polygon is one in which no pair of edges cross each other. A complex polygon 
(sometimes called a non-simple polygon) may have crossing edges. Figure 3(a) gives an 


example of a simple polygon, and Figure 3(b) gives an example of a complex polygon. 


SD Bs 


Figure 3 Examples of simple and complex polygon 


1.4 Current domain 


Filling of complex polygons has been chosen as the application domain, since it has a 
graphically intuitive output. The focus is particularly on raster devices such as printers or 
raster displays, in which pixels within polygonal areas are to be colored. It should be 
noted, however, that the preprocessing of complex polygons into simple polygons 


developed here, could be applied to any of the other general problems mentioned above. 


1.5__Insideness of a point 


Filling a polygon involves finding the points that lie inside the polygon. While the 
question of which points lie inside a polygon is intuitively obvious in the case of simple 
polygons, the question is less obvious in the case of complex polygons. There are two 
rules that are defined to find whether a point lies inside a polygon: the even-odd parity 
rule[13] and the non-zero winding number (NZW) rule[13]. The even-odd parity (EO) 
rule is simpler to implement, but does not intuitively lend itself to self-overlapping areas. 
The NZW rule provides a result, which is consistent with the notion that self-overlapping 


areas should be filled. 


1.5.1 Non-zero winding number rule 


The non-zero winding number rule determines if a point is inside the polygon by walking 
along the entire polygonal path, and counting the number of times the point is encircled 
in a counterclockwise direction (with clockwise encirclements counted as negative). If 
the resulting winding number is non-zero, the point is “inside” the polygon by this rule. 
Note that the polygon may be regarded as a directed graph for this purpose, although the 
orientation is immaterial for the insideness test since all that is required is that the 
winding number be non-zero. 

An alternative way of calculating the winding number is by drawing a ray from 
the point in question to infinity in any direction and then examining all the polygon edges 


that intersect the ray. All polygon edges that cross the ray in a counterclockwise direction 


(relative to the point) contribute +1 to the winding number, and all edges that cross in a 
clockwise direction contribute —1. 

Figure 4(a) shows a complex polygon, star-shaped with crossing edges that was 
filled using this rule. The pentagonal area in the middle lies inside the polygon according 
to this rule, being circumnavigated once, and hence it is filled. The two doughnut shaped 
rings, one with an inner circle with counter-clockwise orientation and the other with a 
clockwise orientation are filled differently, using this rule. A point P is considered inside 
the “polygon” in Figure 4(d). A ray is drawn from it, extended in any direction, intersects 
the polygon at two places, i; and iz. The winding number associated with the point P is 
initialized to zero. When the ray crosses the point 11, the winding number is incremented 
because the direction of the edge it intersected is counterclockwise. Similarly, at iz the 
counter is incremented. The final value of the winding number is non-zero and hence the 
point P is filled. All the other points, enclosed within the inner-circle are filled likewise, 


as shown in Figure 4(d). 


(a) (b) 


(d) 


(c) 


Figure 4 Non-zero winding number rule 


1.5.2 Even-odd parity rule 


The even-odd parity (EO) rule determines whether a point lies inside or outside a 
polygon. This rule does not regard the polygon as a directed graph as in the case of the 
NZW tule, but rather as an undirected graph. A ray is drawn from the point in question in 
any direction and the number of times the ray is intersected by the polygon edges is 


counted. If that number is even the point is outside, otherwise it is inside. 


(a) (b) 


Figure 5 Even-odd parity rule 


Figure 5 shows the filling of the polygons using the EO rule. A point P is 
considered inside pentagonal region of the star-shaped polygon in Figure 5(b). A ray 
drawn from P intersects the polygon at two points, i; and ip which is even. Therefore, the 
region containing P is not filled according to this EO rule. Similarly, the region inside the 
inner circle in Figure 5(c) is considered to be outside using this method. The drawback in 
this method is that, it does not render self-intersecting areas in an intuitive way. For 
instance, the sausage-shaped area in Figure 5(a) created using Microsoft Word does not 
have its self-overlapping area filled because Microsoft Word uses the EO rule to fill its 
polygon. Here, the NZW rule gives a more intuitive result as seen in Figure 4(b). 

Clearly, these two rules differ in their results, and their choice is application 
specific. Most drawing applications use the EO rule because it is simple. On the other 
hand, Postscript, has an operation to select which of the two rules to use, but uses the 


NZW rule as its default. 


1.6 Methods for filling polygons 


1.6.1 Image versus object precision 


In order to characterize current methods for polygon filling, it is useful to introduce the 
notions of image-precision versus object-precision algorithms. An underlying construct in 
image-precision algorithms is that of pixels or scan lines (pixels organized in horizontal 
rows). That is, the algorithm traversal is done on a pixel by pixel, or scan line by scan 
line, basis. The work required of the algorithm is therefore dependent (typically in an 
O(n’) fashion, where n is the number of pixels per unit distance) on the resolution of the 
display raster. On the other hand, object-precision algorithms perform their computations 
in a mathematical manner, independent of any display characteristics. Image-precision 
algorithms are often simpler to design and implement. They use integer arithmetic and 
they trade-off accuracy and speed. In high-resolution environments they provide more 
accuracy, but lower speeds. As device resolutions increase, the speed of these algorithms, 
which may have been adequate for lower resolutions, becomes more of a bottleneck. 
Object-precision algorithms are resolution-independent, but require greater use of 


floating point numbers and hence tend to be slower. 


1.6.2 Current method 


The current technique used in polygon filling is taken from the printer industry (e.g., 
Minolta-QMS), which uses page description languages (PDL), such as PostScript°, HP 
PCL, or HP-GL. A public domain source of sophisticated code for PostScript is 
ghostscript®. These PDLs use, as their method for rendering filled polygons, an image- 
precision algorithm, the scan line algorithm (SLA) adapted for EO or NZW filling, and 
further adapted to generate a list of trapezoids (or triangles, etc.). The SLA algorithm is 


described in Section 1.7. 


1.6.3 Proposed method 


The proposed method for generating a covering set of trapezoids uses two object- 


precision algorithms in tandem; first the current algorithm is used as a preprocessing step 


to convert complex polygons into a set of simple polygons (the left-most arrow in Figure 
1), and this is followed by an already developed trapezoidation algorithm which takes the 
simple polygons generated by the current algorithm and generates the trapezoids (the 


second arrow in Figure 1). 


1.7__Filling polygons using scan conversions 


The most common method of filling polygons in raster devices, in which pixels are 
organized in horizontal rows called scan lines, is by the well-known scan line algorithm. 
A scan line scans every pixel in that line, checks for intersections with the edges and fills 
every pixel, appropriately, by finding whether it is inside or outside the polygon, using 
one of the above-mentioned rules. This process of discovering the pixels that must be set 
is called scan conversion. 

An outline of the SLA follows. Initially a global edge table (GET) is built, which 
contains all the edges of the polygon sorted by their smaller y- coordinate using bucket 
sort, every bucket representing a scan line. Within each bucket, edges are arranged in an 
increasing order of x-coordinate, of the lower endpoint of an edge. Also, every edge 
contains details about the maximum y-value, x-coordinate of the bottom end point and the 
x-increment used to step from one scan line to another. 

Another main data structure is the active edge table (AET), that contains the 
details about the current scan line and the edges associated with this scan line, to facilitate 
further processing using the stored information. 

Thus, for every scan line: 

1. Find all the intersection points of the scan line with the edges of the polygon. 
2. Sort the intersections by increasing x-coordinate. 
3. Fill in all the pixels between pairs of intersections (also called scan run) that lie 


interior to the polygon, using the EO or NZW rule. 


Figure 6 shows scan line filling for a simple polygon. The horizontal scan line at 6, 
intersects the polygon at points 4, 8, 12 and 14 and that is the sorted list of the 


intersection points. 


Scan line at 6 
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Figure 6 Scan Line Algorithm 


1.7.1 Edge coherence 


Calculating the intersections needs to be efficient, so checking every edge of the polygon 
for intersection should be avoided. It is observed that an edge that is intersected by the 
scan line 7, is also intersected by the scan line i+1. This is called edge coherence or scan 
coherence. 

The next x intersection for the new scan line is calculated by using, 


Xie1=xi+1/m 


where m is the slope of an edge. 


1.7.2 Other issues 


There are many other issues involved in filling a polygon and making the result appealing 
to the viewer. For example, an edge that is shared by two polygons, if colored by both the 
polygons, leaves a bizarre color as a result. Hence, it is conventional that, the leftmost 
edge and the horizontal edge of a polygon always belong to the current polygon being 


processed. 


Chapter 2 
RELATED WORK 


2.1 Introduction 


The following is a review of papers that relate to the current research. It includes 
references to research on (simple and complex) polygon partitioning, and line intersection 


problems for complex polygons, or collections of line segments. 


2.2 Partitioning algorithms 


The problem of partitioning a polygon is dated very early in the literature around 1911, 
when the first paper on triangulation of a polyhedron and related problems were 
addressed by Lennes [24]. Several algorithms for polygon partitioning in the forms of 
triangulation, trapezoidation have been reported in the literature since 1978. Monotone 
partitioning by Garey, Johnson, Preparata, and Tarjan (1978) [16] is an algorithm that 
partitions a simple polygon into monotone polygons, and extends the concept to carry out 
triangulation. Since then, there have been numerous algorithms devised and implemented 
that sought solutions for triangulation, with less run-time complexity [14][28][11]. 
Finally, in 1991, Chazelle [9] came up with an algorithm that triangulates a simple 
polygon in linear time. But this algorithm was claimed to be ‘unimplementable’, because 
of its enormous details. All the above-mentioned algorithms were written to partition a 
simple polygon into a set of triangles, using various techniques. 

In 1995, Bruce Romney [27] used the concept introduced by Chazelle and Incerpi 
to simplify the required data structures, to vertically decompose the faces of an 


arrangement into trapezoids, which was extended for triangulation. Also, there is Seidel’s 


10 


algorithm for triangulation, which is simpler, but not a faster triangulation algorithm that 
can be widely used in the industry. Hain’s algorithm [17] for trapezoidation is a more 
practical solution for partitioning, which is almost linear in complexity except in certain 
conditions. 

More recently, Nancy M. Amoto, Michael Goodrich and Edgar Ramos [2] 
published an algorithm for triangulation of a simple polygon in linear time, through 
trapezoidation of a polygon. They have also stated that it is simpler than Chazelle’s 
optimal deterministic algorithm. 

All these algorithms have varied application domains. As mentioned earlier, the 
domain of interest for this thesis is filling the polygons. This application can be used in 
the printer industry and various graphics applications, where even some complex 
characters are considered as polygons to be filled. In this context, problems involved like 
point location, etc., are discussed by Lee and Preparata in their papers on the applications 
involved in point location on a planar subdivision [23][22] in 1977, before the break- 


throughs for triangulation problems that were based on these ideas. 


2.3_Line intersections 


The other main issue involving polygon partitioning, i.e., finding the intersections 
between collections of line segments, was addressed by Bentley and Ottman in 1979 [6]. 
It was further discussed by Chazelle and Guibas [10] in 1989, which is one of the most 
commonly referenced papers in algorithms that deal with polygon triangulation or 
trapezoidation. A recent paper by Chazelle in 1992 [7], discusses the solutions putforth to 
date, for the line intersection problem and the issues involved in that problem that 
complicate the solutions. It is one of the main referenced papers of this thesis, as it gives 
an overall idea of all the methods involved in finding the solution, for practically all the 


test cases, that might occur in the application domain. 


2.4 Other sources 


The above-mentioned papers do not represent any particular area for application. The 
Postscript manual documentation [1] describes the process of filling the polygons that is 


followed as standard in the Postscript specification. Also the standard methods, described 
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in the books on primitive methods in computer graphics on rasterization, were read and 
noted [13]. 

After a careful search in the literature, it was found that no algorithm has been 
written to partition a non-simple polygon into a set of simple polygons. Most of the 
algorithms, just avoided the complex polygons, and processed only the simple polygons. 
This research proposes to fill-in the gap, by providing an implementation for the two 
main algorithms outlined in Chapter 4. 

The algorithms used for filling polygons that have been commonly used by printer 
companies (like Minolta-QMS), and Page Description Languages, like PostScript [1] 
were considered for a comparison study in this thesis, but the code for those was not 
available. However, it is known that the graphics library uses the SLA algorithm for 
polygon partitioning [18]. 

In particular, the code for a sophisticated public domain version of a PostScript 
interpreter and renderer, ghostscript [15], was studied and found to also implement 


polygon partitioning by the SLA algorithm. 
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Chapter 3 
METHODOLOGY 


3.1 Aspects of work 


This work concerns a new algorithm to partition a non-simple (complex) polygon into a 
set of simple polygons, using either the NZW rule or the EO rule, thereby providing a 
preprocessor for many applications involving polygon partitioning. The work includes 
the following components 

1. Design the algorithm based on an outline provided by Dr. Hain. 

2. Implement the algorithm. 

3. Create a tool to interactively generate and manipulate polygons (simple and 

complex). This is useful for the dynamic testing of the algorithm. 
4. Test and demonstrate the correctness of the algorithm using this tool. 
5. Provide a test bed to demonstrate the relative performance of the algorithm 


compared to current methods. 


The algorithm clearly should be correct, and should be efficient within useful, i.e., 
practical domains. The algorithm will be empirically compared to the SLA algorithm, the 


algorithm currently used in industry. 


Implementing the algorithm involves determining the data structures that would 
best solve the problem, minimizing performance overhead. Some ideas found in the 


literature were adapted for this thesis. 
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3.2 Correctness 


Correctness testing required the generation of input polygons, and availability of their 
correctly filled outputs. Two generation methods were used in this work, and embodied 
in a single interactive tool built for this purpose: 

1. Random polygon generation. 

2. Interactive polygon generation. This method of testing is found to be useful in 
testing a variety of computational geometry algorithms because of their visual 
outputs. Because of the dynamic nature of this type of testing, a very large domain 
of equivalence classes can be covered. This technique was described in [18]. 
Code was adapted from [4] [5] to generate random polygons for a given number 

of vertices. This random polygon generator (RPG) is an application developed as a part 
of Thomas Auer’s Master’s thesis. It is a set of C programs that use an X-Windows API. 
This application has few heuristics, and can be utilized to generate different kinds of 
polygons, like star-shaped, monotone, etc. These polygons are random in their geometric 
properties, with no restriction on their sinuosity (winding property,) or their horizontal or 
vertical co-linearity. 

Interactively produced polygons can be generated and edited using the interactive 
testing tool. The tool allows the user to move, add, or delete vertices, starting with either 
a simple triangle, or a randomly generated polygon. The user can generate multiple 
polygons, within the drawing area. 

These polygons can be saved as either a text file, or as a PostScript program. The 
PostScript output is viewable with a PostScript viewer (ghostscript)—which is presumed 
to provide the correct output—to enable visual comparison with the results from the 
NZW algorithm. The correctness was taken as a consistency between the results from the 
PostScript viewer, the SLA and the NZW algorithms, together with, in relatively simple 
cases, manual solutions. 

The tool showing a sample complex polygon with a randomly generated 20- 


vertices polygon, is shown in Figure 7. 
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Figure 7 Interactive testing tool showing random input polygon 


3.3. Performance 


For performance testing both the algorithms (NZW and SLA) were embedded into a 
Windows console application, and used a high resolution (sub-millisecond) timer 
provided within the MS Windows® operating system. The timing was performed on a 
release compile (i.e., optimized, and without debug code.) All possible background 
processes were removed to make sure that as many of the resources of the computer were 
devoted to the operation of the algorithms. The timing was restricted to only those 
operations performed by the algorithms necessary to determine the outline of the fill 
operation. No scan line spans were actually filled in SLA. No actual rendering was 
performed. Care was taken to avoid including any IO overhead within the timed brackets. 
The timing data suite was a set of randomly generated complex polygons 
classified by their number of vertices, number of intersections, and height in pixels. The 


generation of this data suite is further discussed in Section 5.2.2.1. 
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3.4 Hypotheses 


1. An object-precision algorithm can be devised to partition a complex polygon into a 
set of simple polygons using the NZW rule (NZW Partitioning algorithm.) 

2. The devised algorithm will generate a filled polygon no slower than, and, in more 
complicated cases, faster than the SLA. 

Hypothesis | can be verified by testing the correctness of the algorithm, as outlined in 

Section 3.2. Hypothesis 2 can be verified by comparative performance measurements, as 


outlined in Section 3.3. 
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Chapter 4 
THEORY 


In this chapter we describe a new algorithm to partition a complex polygon into a set of 
simple polygons which outline the filled area. This partitioning is with respect to a 
particular fill rule—either the non-zero winding number rule described in Section1.5.1, or 
the even-odd parity rule described in Section1.5.2. 

The initial investigations led to the conclusion that the EO partitioning algorithm 
could be easily designed and implemented using the primitives developed to design and 
implement the NZW partitioning algorithm. Thus, the focus was to outline an algorithm 
to partition a complex polygon into a set of simple polygons (possibly with holes) using 
the NZW rule. 

At present, the algorithm presented does not address any singularities such as 
degenerate polygons (e.g., zero size, coincident vertices, etc.), with the assumption that 
the algorithm can be extended to include these cases without loss of generality, if the 


fundamental concept is correct. 


4.1 Outline of NZW partitioning algorithm 
We will first outline the NZW partitioning algorithm. It should be noted that 
many implementation details will be omitted at this point, and only a high-level 
understanding of the algorithm will be provided in this section. We will use the notion 
that the interior of a simple polygon can be considered as a (an infinite) set of points in a 


two-dimensional space. 
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4.1.1 Terminology 


This section formally defines the terminology used in the description of the 
outline of the NZW algorithm. 

The input to the algorithm is an ordered list of vertices—that we denote as 
polygon vertices—of a (complex) polygon. The ordering of the polygon vertices induces 
directed line segments called polygon edges. 

An intersection is defined as the point at which exactly two polygon edges 
coincide. When more than two polygon edges intersect at the same point, multiple 
coincident intersections result since each intersection results from a pair of edges. 

A directed polygon edge may be divided at intersection points into two or more 
collinear, directed (in the same sense as the parent edge) edge segments. Each polygon 
edge segment may be thought of as belonging or being the child of a polygon edge. Both 
the in- degree and out-degree of an intersection is two. The described algorithm considers 
the intersection as having two coincident pseudo-vertices which are companion to each 
other. For the sake of orthogonality, we will consider a polygon vertex as an intersection, 
and define two associated pseudo-vertices, one of them being null. When we wish to be 
specific, we will distinguish the pseudo-vertices as polygon pseudo-vertices, and 


intersection pseudo-vertices. 


infinitesimal 
gap 
v 
xX Se XI X2 
Edge 
Edges segments 
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Figure 8 Pseudo-vertices x1, and x2 at intersection point X 
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4.1.2 Algorithm overview 


The algorithm joins, at each pseudo-vertex, an incoming edge segment belonging to one 
polygon edge with the outgoing edge segment belonging to the other polygon edge, as 
shown in Figure 8. By connecting the pseudo-vertices in this way, a set of polygons are 
formed which we will term as component polygons. Thus, the edges of component 
polygons are the edge segments of the complex polygon, and the vertices of the 
component polygon are the pseudo-vertices. Every pseudo-vertex will be associated with 
a component polygon, and every component polygon will contain a disjoint set of 
pseudo-vertices. The following theorems characterize component polygons. 
Theorem 1: Component polygons are simple. 
Proof: No edge segments intersect, otherwise they would be partitioned at that 
intersection. Each pseudo-vertex is traversed exactly once. If we conceptually 
consider an infinitesimal gap between the companion pseudo-vertices, as in 
Figure 8, it can be seen that the sequence of directed edges does not cross any gap 
between pseudo vertices, and therefore cannot pass the same intersection more 
than once, and therefore cannot self-intersect. Since there is a finite number of 
pseudo-vertices, the sequence must close because of the pigeon hole principle, 


and therefore forms a simple polygon. a 


Theorem 2: Component polygons are disjoint, or fully enclose, or are fully 
enclosed by another component polygon. That is, in a set of component polygons 


P, 1<is<n, for each pair P and P, withi # j one of the following is true: 


le aes, 

2. BoP, 

3. PAP, =O 
Proof: 


It is seen that no component polygon can intersect with another component 
polygon at either a component polygon edge (edge segment) or, again considering 
a conceptual gap between pseudo-vertices, at a component polygon vertex 


(pseudo-vertex). Thus the component polygons can only be nested or disjoint. um 
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Figure 9 (a) complex polygon with labeled polygon vertices 
(b) the result of the filling by NZW algorithm 


Each component polygon, being simple, is oriented either counter-clockwise or 
clockwise, and so contributes a winding number of +1, or —1 to any enclosed point. We 
therefore associate a winding number of +1 with each component polygon (or 
alternatively with each of its pseudo-vertices.) To determine whether a given point is 
“inside” the original complex polygon, one may simply add the winding numbers of all 
the component polygons that enclose the point, to see if it is non-zero. We will call the 
sum of winding numbers of the enclosing polygons the net winding number. Since the 
component polygons are nested, we can associate a net winding number with each 
polygon (or alternatively with each of its pseudo-vertices.) We can also visualize the 
component polygons as stacking together with the smaller polygons “on top of” the 
largest one. The net winding number of any point is the sum of the winding numbers of 
all the polygons stacked at that point. 

We will define a vertex of a simple polygon as being either convex (right-turning 
for a clockwise oriented polygon, or left-turning for a counter-clockwise polygon), or 


concave (left-turning for a clockwise polygon, or right-turning for a counter-clockwise 


polygon). 
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The algorithm proceeds by finding the left-most polygon vertex, and 
corresponding polygon pseudo-vertex. This pseudo-vertex will necessarily be convex, 
and will define the orientation of the associated component polygon, and consequently 
the net winding number of this pseudo-vertex can be initialized to either —1 or +1. 

As the component polygon is traversed starting at the identified pseudo-vertex, 
every companion to a convex pseudo-vertex is added to a queue (alternatively, a stack 
can also be used). This companion will be the starting point for the traversal of a 
(potentially) new component polygon disjoint to the current one. Actually, this disjoint 
polygon may “touch” the current polygon at several places. Since we know that the 
disjoint polygon will have the opposite orientation as the parent, we can update the net 
winding number of the companion appropriately. 

Similarly, as the traversal proceeds, all companions to a concave pseudo-vertex 
will be added to the queue, and will be the starting point for the traversal of a 
(potentially) new component polygon that will be an outermost component nested within 
the current one. Again, this inner component polygon may touch the current component 
polygon at several places. Since we know that the disjoint polygon will have the same 
orientation as the parent, we can update the net winding number of the companion 
appropriately. 

Once the current component is completely traversed, a pseudo-vertex is dequeued, 
and, if it has not already been traversed, starts a new traversal in a like manner to the 
previous one. This process continues until the queue is empty. 

At this point we have a collection of all the component polygons, each with an 
accumulated net winding number. The outline of the area that is to be filled is the subset 
of component polygons having net winding number of —1 or +1 (forming the outer 
boundaries of the filled areas), and component polygons having a net winding number of 


0 (forming the inner boundaries, or “holes”.) 


4.13 Example 


Let us consider the input complex polygon in Figure 9. This example figure has most of 


the special cases that need to be handled by the algorithm. Also, it should be noted that 
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the input polygon edges are directed, which is a factor that decides the insideness of a 
point in a polygon in the NZW filling rule as seen in Section1.5.1. 

The component polygons of the example polygon in Figure 9 are shown in Figure 
10. The first component polygon generated in the algorithm is the outer-most polygon P 
(Figure 10(a)), which has a winding number of -/ , filled with downward slanting lines, 
as shown. The second and third polygons P» (Figure 10(b)) and P3 (Figure 10(c)) have a 
winding number of —1 in accordance with their orientation, and have net winding 
numbers of —2 (since P2 overlays P,) and -3 (since P3 overlays P; and P2) respectively. 
The fourth polygon P,4 in Figure 10(d) is not filled because its net winding number is 
Zero, since it has an opposite orientation to its encompassing polygon, P; and hence it 
forms a hole inside P;. Thus, holes are handled in a natural way by the NZW algorithm. 
The fifth polygon Ps in Figure 10(e) is disjoint from P; and has a winding number of +1. 
Figure 10(f) shows the overlapping component polygons, and each area with different 
shading has a different net winding number. 

When the input complex polygon is processed by the NZW algorithm, three 
simple polygons define the filled area as shown in the Figure 9(b), and one of these is a 


hole. 
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Figure 10 (a) -(e) Component polygons of the example input polygon. (f) 
Overlapping component polygons 


4.2 Implementation Details 


4.2.1 Intersection points 


Each intersection between a pair of polygon edges is calculated analytically in this 
algorithm using parametric equations. Naturally, each polygon edge may intersect with 
up to n-1 other edges, and we want to associate this list of intersections with either that 
edge or, say, its starting vertex. 

The evaluation of all intersections is O (n*) in the worst case, since every pair of 
polygon edges potentially intersects. The related problem of finding the intersection 
points of a collection of line segments was addressed by Chazelle [7], and is itself a 
surprisingly difficult task to do optimally, and we do not propose to implement a 
theoretically optimal algorithm. Nonetheless, the discovery of intersections probably 
represents the major computational work of the algorithm proposed here. However, with 
reasonable heuristics, and an assumption that input polygons are not extremely 


convoluted in practice, the computation can be considerably reduced. We tried to 
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minimize the need for intersection calculations by taking advantage of any occurrence of 
monotonic chains (in either the x or the y direction), but the overhead was found to be 
greater than the gains in most typical situations. Trivial rejection of line segment 
intersections by non-overlapping bounding boxes was used to avoid intersection 


calculations, and was found to be advantageous. 


4.2.2 Intersection Linkage 


When the input polygon is traversed to generate the component polygons, each pseudo- 
vertex that is traversed, has a pointer to the next pseudo-vertex it is linked with. This 
information that is contained in every pseudo-vertex is called its intersection linkage. The 
following paragraphs briefly introduce the classes that were constructed to establish the 
necessary linkage. It should be noted that the objects created with these classes, are 
temporary, i.e. they are needed until the intersection linkage is complete. 

A polygon edge array of polygon edges with each edge pointing to a list of 
intersections, sorted by their parametric value relative to the edge (i.e., t = 0 at the start 
vertex, t= 1 at the end vertex of the edge). This task will give rise to the companion 
pseudo-vertices, one for each intersecting line segment. The polygon edge array that will 
be generated for the example in Figure 9 is shown in Figure 12. 

An intersection object is composed of a pair of pseudo-vertices at an intersection. 
This includes both the polygon pseudo-vertices and the intersection pseudo-vertices. The 
polygon edge array created above is used to create an intersection master list that 
contains all the intersection objects. The intersection master list that will be generated for 
the example in Figure 9 is shown in Figure 16(a). The polygon edge array and the 
intersection master list are built simultaneously and their data structures are discussed in 


Section 4.2.3.2.1 


4.2.3. Data Structures 


The challenge in the implementation of the algorithm lay in the design of the appropriate 
data structures. In this section, the data structures of all the objects required to implement 


the algorithm are illustrated in Figure 11 and Figure 13. 
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4.2.3.1 Support classes 

The basic entities that are fundamental are as follows: A polygon vertex is stored as an 
object of type Vertex, aclass containing two float values, for x and y co-ordinates. 
The input polygon is stored as an object of type Polygon, a class containing a set of 
ordered vertices, implemented as a linked list of vertices, a vertex list. A polygon edge is 
stored as an object of type LineSegment, a class containing a start vertex and an end 
vertex. With all these supporting classes MultiPolygon (more than one polygon e.g., a 
set of component polygons, or a set of outline polygons together with hole polygons) is a 
class that contains a linked list of polygons ( poly list ), with each polygon in turn 
containing a pointer to a vertex list. This is illustrated in Figure 11. The linked list, vector 


and queue were implemented using the container classes from STL. 


Vertex LineSegment 
float x, Vertex start; 
floaty; Vertex end; 

(a) (b) 

Polygon MultiPolygon 

list<Vertex> list<Polygon> 
vertexList; polyList; 

(c) (d) 


Figure 11 (a) Vertex object (b) LineSegment object, (c) Polygon object and (d) 
MultiPolygon object 
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Figure 12 The Polygon Edge Array for the Figure 9(a) 


4.2.3.2 Intermediate classes 


In order to get the master list of intersection objects that contains the necessary 
information, various intermediate objects were designed. The polygon edge array was 
implemented using an object of a class called nVertex which contains a vertex 
(polygon/intersection pseudo-vertex), the intersecting line segment(/), the parametric 
value (t¢ ), with respect to the line segment and the index of the position of the pseudo- 
vertex, in the intersection master list, as shown in Figure 13(b). The array of polygon 
edges is implemented as a vector (nVertex objects holding polygon vertices), each having 
a pointer to a linked list of intersections (nVertex objects of intersections), as shown in 
Figure 12 . It is implemented with the nVertex objects that hold necessary information to 
create the master list. 

The most important object is an intersection object, which contains information about 
every pseudo-vertex in the input polygon. Every intersection object contains an 
intersection point (v), the starting vertices of the two intersecting polygon edges (origin 
vertex! and origin vertex2), the pointers to the next intersection objects i.e., the value of 


the index position in the intersection master list (index/ and index2) and the winding 
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number of the polygon generated by the corresponding polygon vertex (winding number), 


as shown in Figure 13(a). 


nvertex Intersection 
Vertex Vv; Vertex Vv; 
float parameter; Vertex originvertex 1; 
int index; Vertex originvertex2; 
Line Segment /; int index); 
int index»; 
int windingnumber; 
(a) (b) 


Figure 13 (a) nVertex object (b) Intersection object 


4.2.3.2.1 Setting up the intersection linkage: 

The intersection master list and the polygon edge array are built simultaneously, 
since the intersection master list is required to fill the index field for the objects in the 
polygon edge array. Once all the indices are filled in the polygon edge array, then the 
values in the intersection master list are filled from the information available from the 
polygon edge array. For every intersection object in the intersection master list: 

1. search for the vertex in the polygon edge array. 
2. update the index; and index? value of the intersection object, from the 
corresponding object in the polygon edge array. 

If the vertex is a polygon vertex, then it is searched in the polygon edge array. 
There is only one index to be updated in the intersection object. The other index is 
marked -1, to imply that the path is invalid. If it is an intersection, then it is searched in 
the intersection list of every polygon edge, in the polygon edge array. Since an 
intersection object also has information about intersecting polygon edges, the index, 
corresponds to origin vertex; and index2 corresponds to origin vertex. 

For illustration , let us consider the example polygon in Figure 9. The intersection 


object containing the pseudo-vertex 7; in the master list will have the data as shown in the 


Zi 


Figure 14(b). The pseudo-intersection vertex i, is linked as an intersection from both the 
edges p, and pz, in the polygon edge array, as illustrated in the Figure 14(a). When it is 
linked under the edge p,, the index position of i; denotes the next intersection object in 
the intersection master list it is linked with, with respect to the corresponding polygon 
edge, which is i7. Similarly, with the polygon edge starting at pz, it is linked to the 


intersection which happens to be the polygon vertex ps. 
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Figure 14 (a) is the illustration of i; in the polygon edge array (b) is the intersection 
object of i; in the intersection master list 


The data structure for a multi-polygon has also been discussed, because the 
algorithm can successfully handle more than one polygon, as long as there is at least one 
intersection point between them. In this case, the intersection points are also calculated, 


between the polygon edges of every polygon with each other. 


28 


4.2.3.3 Illustration of the intermediate classes using an example 
This section consists of a more detailed view of the algorithm by illustrating the objects 
that are created from the intermediate classes. For this purpose the same example in 


Figure 9 is considered. 


After calculating all the intersection points between the polygon edges (of one or more 
polygons), the edge segments are generated. The intersection points and the edge 
segments are shown in Figure 15. All of the intersection points and the polygon vertices 


are then converted into pseudo-vertices. 


Figure 15 Complex polygon showing edge segments and pseudo-vertices 


A polygon edge array, with a list of intersection vertices, sorted according to its 
parametric value, ¢, is created as shown Figure 12. The intersection master list that has the 
intersection objects with appropriate pointers, required to traverse the polygon is 
generated. The list generated for this particular example is shown in Figure 16(a). For this 
example, the starting pseudo-vertex with the lowest co-ordinate tuple is the polygon 
vertex p;. Meanwhile the intersection objects that are visited are pushed into the 
intersection queue, which is shown in Figure 16(b) for the cycle generated. The indices of 
the objects that participated in this cycle are updated (exhausted) in the intersection 


master list. 
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The next starting pseudo-vertex is identified by searching the intersection queue. 
For this example it is at i79 as shown in Figure 16(b). The complete cycle that generates the 
next nested component polygon with the vertices. This procedure is repeated until the 
intersection queue is empty. The values in the intersection queue during the traversal are 


shown in Figure 16(b), for few component polygons. 


P1| P2| P3 | P4| Ps | Po| P7 | Ps | Po | Pro | t1 | i2 | 13 | t4 | is | t6 | 7 | ts | to | t10 


11 | 12 | ig | is | ig | iro | i7 | ig | io | iy | pr | iy | ps | iz | iz | ig | po | is | pa | tro | Pro 


(b) 


Figure 16 (a) intersection master list of pseudo-vertices (b) Shows various values in the 
intersection queue 
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Chapter 5 
RESULTS 


The NZW algorithm for the partitioning of a non-simple polygon into a set of simple 
polygons was implemented using the data structures explained and illustrated in section 
4.2.3. An implementation for the NZW algorithm was built using the C++ language, 
using a class collection provided by Dr. Hain. An interactive testing tool to view the 
output of the algorithm was implemented using C++ and MFC (Microsoft Foundation 
Classes, and object-oriented set of classes encapsulating the MS Windows API.) Finally, 
a timing framework was built using a Windows console application, and a high resolution 


timer. 


5.1 Correctness 


Correctness testing required the generation of input polygons. Two methods were used in 
this work, and are described in detail in Section 3.2: 

1. Random polygon generation. 

2. Interactive polygon generation. 

The algorithm showed consistency between the SLA algorithm, the NZW 
algorithm, and selected polygons that were compared with the ghostscript output. The 
only precondition for the algorithm is that the polygons are not degenerate (zero area). It 
should also be noted that the algorithm works for many complex polygons, provided 
there is at least one intersection point between them. Hence the robustness of the 
algorithm is demonstrated. Note the difference between the results of filling using the 


NZW and the EO rule, as shown in Figure 17, Figure 18 and Figure 19. 
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Figure 17 NZW fill using NZW filling rule 


Figure 18 SLA fill using NZW filling rule 
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Figure 19 The SLA fill using SLA filling rule 


5.2 Performance 


The performance of the NZW algorithm was compared with the scan line algorithm, 
which is the current filling technique typically used in industry. For this purpose, the SLA 
was implemented’ based on an outline given in [13]. The following sections dicuss both 


the empirical and the theoretical correctness of the NZW algorithm. 


5.2.1 Theoretical complexity 


The factors affecting the performance of the NZW and SLA algorithm are assumed to be: 
1. the number of vertices (7), 
2. the number of intersections (7), and 
3. the screen resolution (height in pixels) 
The theoretical complexity of the NZW algorithm is O(s logs) where s is the number of 
intersections (including the polygon vertices, also regarded as intersections) in the input 
complex polygon. Since the algorithm basically traverses intersection points (a linear 


contribution,) and sorts intersection on an edge, with the worst-case being all 


' Both the filling rules, NZW as well as EO were provided. 
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intersections occurring on a single edge, yielding O(s logs). The number of intersection 
is bounded by O(n’), where n is the number of vertices (i.e., edges) in the input polygon. 
Thus, in the worst case, the NZW algorithm has complexity O(n’ logn). Since this 
algorithm is object-precision, the height of the polygon in pixels, which is dependent on 
the screen resolution, should not affect its performance. 

On the other hand, the SLA is an image-precision method, which scans every scan 
line from top-to-bottom and left-to-right. So the time taken to fill a complex polygon by 
the SLA, using the non-zero winding number filling rule was assumed to be directly 
proportional to both the width and the height of the polygon. The SLA calculates the 
intersection of the polygon edges with the scan line, while scan-converting, using edge- 
coherence. Thus, the SLA does not partition a polygon, but only decides on which parts 
need to be filled. SLA does not require the intersections points to be stored, but a test of 
whether the intersection of two edges with the current scan line changes their relative 
left-to-right order (i.e., two edges intersect) must be made in any case. Any change in the 
number of intersections, should have no effect on the SLA, for a given number of vertices 
of a polygon and for a given height of the polygon in pixels. 

In addition, we also wanted to study the performance of the algorithms with 
respect to the number of vertices of the polygons, for a given number of intersections and 


height. 


5.2.2. Empirical relative performance 


5.2.2.1 Performance test suite 
A test data set was constructed that consisted of 6 groups of randomly generated 
polygons whose vertex count were 20, 30, 40, 60, 80 and 100, and each of those having 4 
different numbers of intersections (depending on the vertex count). These 24 polygons 
were again categorized into six separate groups, scaled with different heights of (in the 
order of) 200, 400, 800, 1600, 3200 and 6400 pixels. 

The base polygons were first generated for the given number of vertices and 


intersections, in a rectangular area of 420 x 700, with a height of 400 pixels. From these 
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base polygons, the required data sets were then generated by scaling the y-values by 0.5, 
2,4, 8 and 16, respectively. The scaling was done only in the y-direction, because the 
SLA depends on the number of scan lines and not in the x-direction’. There are 10 
distinct polygons in each category, and this is believed to provide a sufficiently 
representative sampling. 

The different number of vertices, number of intersections and the heights of the 
polygons were selected based on the discussion with Dr. Hain, who is quite 
knowledgable in this field. It was gathered that most commonly used polygons for 
printing had vertices in the range of 40-100, with less than 50 intersections and their 
average height was 6400 pixels. The height chosen is a representation of any figure that 
around 5” long (at 600 or 1200 dots per inch, equivalent to respectively 3200 and 6400 
pixels) on a 1200 dots per inch high-resolution printer, which is a reasonable dimension 
for comparison. Since we do not have an actual data set of practical polygons, we 
compared the algorithms with polygons of typical size that are used in the printer 
industry. 

For comparison, both the algorithms, NZW and SLA were executed with the same 


set of data on a machine with the system specifications as shown in Table 1. 


Table 1 System Specifications 


OS Name Microsoft Windows XP Professional 

Version 5.1.2600 Service Pack 1 Build 2600 

OS Manufacturer Microsoft Corporation 

System Manufacturer Gateway 

System Model E-3600 

System Type X86-based PC 

Processor x86 Family 15 Model 1 Stepping 2 
GenuinelIntel ~1694 Mhz 

BIOS Version/Date Intel Corp. 
PT84510A.15A.0003.P01.0110311619, 


* No actual filling was done, and so the algorithm is independent of the number of pixels in the x-direction. 


35 


10/31/2001 

SMBIOS Version 23 

Hardware Abstraction Layer Version = "5.1.2600.1106 (xpsp1.020828- 
1920)" 

Total Physical Memory 512.00 MB 

Available Physical Memory 251.50 MB 

Total Virtual Memory 1.67 GB 

Available Virtual Memory 1.10 GB 

Page File Space 1.17 GB 


5.2.2.2 Performance results 

The execution times for the compiler-optimized SLA and NZW algorithms operating on a 
variety of inputs were collected. Specifically, the effect of the number of vertices, the 
number of intersections, and the polygon height on the execution time of both the 


algorithms was measured. Each of these measurements will be discussed separately. 


5.2.2.2.1 Performance versus polygon height 

The effect of height on the execution time for partitioning a polygon with given a 
number of vertices and intersections is shown in Figure 20. It can be seen that for NZW 
the running time is essentially constant with height, as was expected since this is an 
object-precision algorithm. The results of the exection times for three different sets of 
polygons are shown to compare the results with the practical sizes. It can be seen that for 
polygons with 106 vertices and 100 intersections, which are considerably complex, 


NZW has an impressive execution time as compared to SLA. 
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Error! Not a valid link. 


Figure 20 NZW time vs. Height of the polygons (trend lines) 


Conversely, the SLA timings increase linearly with height, again as was expected, 
since the number of operations is proportional to the number of scan lines (i.e., height). 
The graph shown in Figure 21 shows the trend lines of two different sets of polygons, 
with respect to the height of the polygon. It can be seen that, although for polygons of 
height 200 pixels and with around 20 vertices SLA seems to perform as well as NZW, for 
polygons of practical height NZW runs considerably faster than the SLA. For this reason, 
the performance of the NZW in the following sections are discussed for polygons with 


height of 6400 pixels. 


Error! Not a valid link. 


Figure 21 SLA time vs. Height of the polygons (trend lines) 
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5.2.2.2.2 Performance versus number of vertices 


The effect of the number of vertices on the performance of the SLA and NZW 
algorithms was tested from the results obtained. For each of the chosen number of 
vertices, the random polygon generation algorithm was used to produce a set of 10 
polygons. Each of these polygons that were generated had 50 intersections. 

Figure 22 shows the effect of the number of vertices on the execution time for 
polygons, with number of intersections and height kept constant. It can be seen from the 
results that NZW algorithm runs faster than SLA for a sample data of practical size and 
height. 


Error! Not a valid link. 


Figure 22 Performace vs. Number of Vertices (50 intersections) 
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5.2.2.2.3 Performance versus number of intersections 
Timing tests were conducted to find the effect of the number of intersections on 
the NZW and SLA algorithms, keeping the number of vertices constant at 100. The 


results are shown in Figure 23. 


Error! Not a valid link. 


Figure 23 Execution Time vs. Number of Intersections (for 100-vertex polygons) 


It can be seen that the execution time for SLA is relatively independent of the 
number of intersections, as was expected. Surprisingly, this was, to a lesser degree, also 
true for NZW. That is, the NZW algorithm performs well, even with a greater number of 
intersections. Since NZW runs much faster than SLA, it is more efficient than SLA, 


especially for polygons of practical sizes. 
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5.2.2.3 Relative Performance of NZW 

The relative performance (ratio of the execution time of SLA vesus the execution 
time of NZW) as a function of the number of intersections has also been calculated. It is 
shown in Figure 24 for polygons with 20 vertices, and in Figure 25 for polygons with 86 
vertices. Clearly, the greater the number of intersection, the more “complicated” is a 
polygon. It can be seen that SLA runs almost as fast as NZW for quite complicated 
polygons (86 verticecs and 100 intersections) at very low resolution (200 pixels). 
However, for similarly complicated polygons in more demanding scenarios required by 
the higher resolution of modern printers (1200 dot per inch), NZW typically runs more 
than 10 times faster than SLA. 


Error! Not a valid link. 


Figure 24 Performace Ratio of NZW vs. Number of Intersections (20-vertex polygons) 
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Error! Not a valid link. 


Figure 25 Performance Ratio vs. Number of Intersections (86-vertex polygons) 


5.3. Conclusions 


The design and implementation of the NZW algorithm, and the thorough testing of the 
algorithm using the interactive testbed described in Section5.1 confirms the first 
hypothesis which states that an object-precision algorithm can be implemented to 
partition a complex polygon into a set of simple polygons (using the NZW rule.) 

The second hypothesis states that the NZW algorithm fills a complex polygon no 
slower than the SLA, and in complicated cases faster than the SLA. This hypothesis is 
confirmed by the results shown above, and for the following cases: 

1. With the height and number of intersections of polygons kept constant, NZW 

executes faster than SLA for polygons with varying number of vertices, as 


shown in Figure 22. 


4] 


2. With height and number of vertices of polygons kept constant, NZW executes 
faster than SLA for polygons with varying number of intersections, as shown 


in Figure 23 


Section 5.2.2.3 on relative performance as a function of the number of intersections, 
mapped for different heights, further illustrates this hypothesis. The results clearly 
indicate that NZW runs 10 times faster than SLA for practical sizes. This is an amazing 
and gratifying result. Also, since NZW is an object-precision algorithm, there is the side- 
effect of higher precision results. The performance evaluation of the NZW algorithm 
described here shows it to be superior to the most-commonly used algorithm(SLA), and 
makes it a valid and valuable pre-processing algorithm useful in graphics libraries of 


modern laser printers. 


5.4 Future work 


Future work leading from this research could include the following two extensions. 

1. Adapting the NZW algorithm to fill complex polygons using the EO rule. This is 
a relatively simple modification. 

2. Adapting the NZW algorithm to generate the union, intersection, and difference 
polygons of two or more given overlapping simple concave polygons. While 
object-precision algorithms for these operation on convex polygons are well- 
known, the NZW algorithm provides the basis for the same operations on concave 
polygons, which are currently performed using a variation of SLA. The algorithm 
is in fact very close to what is needed. For example, the algorithm is already 
capable, given two overlapping simple concave polygons having the same 
orientation, of producing a component (in almost the same sense as is used in this 
thesis) polygon, with possible holes, representing the union of the two polygons. 
By reversing the orientation of one of the polygons, component polygons 
representing the differences are generated, and the union of these is the inverse of 


the intersection 
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APPENDIX 


APPENDIX A 
CODE 


// Multipolyon.h 
// Written by T. Hain. Extended by L. Subramaniam, Spring, 2003 


ifndef _MULTIPOLY 
define _MULTIPOLY 
define — MFC 


include <list> 
include <vector> 
include <queue> 
include <ctime> 
include <algorithm> 
include <stack> 
include "stdafx.h" 


using namespace std; 
typedef float Coord; 


ZITTLLTTLTTTTATT TITTLE 


// Vertex: vertex 


// 

// Constructors 

// Vertex () 

if Vertex(const int x, const int y) 

// Vertex(const CPoint &point) 

// 

// Assignment 

// void operator= (const Vertex &v) 

// 

// Cast to CPoint 

Lf operator CPoint () 

// 

// Tests if vertex is within distance 'resolution' of line (p1,p2) 
// BOOL isNearLine(const Vertex &pl, const Vertex &p2, int resolution) 


class Vertex: public CObject 
{ 
public: 

Coord x, yi 


Vertex(const Coord _x = 0, const Coord _y = 0): x(_x), y(_y) {} 
Vertex(const Vertex &v): x(v.x), yl(v.y) {} 


void operator= (const Vertex &v) 
{ X =v.x; y=v.y; } 


Vertex operator+ (const Vertex &v) const 
{ return Vertex(x + v.x, y + v.y); } 


48 


Vertex operator- (const Vertex &v) const 
return Vertex(x - v.x, y - v.y); } 


Vertex& operatort+= (const Vertex é&v) 
MPS VG 


y += V.Ve 
return *this; 


Vertex& operator-= (const Vertex &v) 
ST Mae ky 
Yow e Meye 


return *this; 


bool operator== ( const Vertex &v) const 
return x == v.x && y == v.y; } 

bool operator!= (const Vertex &v) const 
return x !=v.x || y !=v.y; } 


bool operator< (const Vertex &v) const 
return x < v.x || x ==v.x && y < v.y; } 


bool operator<= (const Vertex &v) const 
return x <= v.x; } 


bool operator> (const Vertex &v) const 
return y < v.y || y ==v.y 7 } 


#ifdef | MFC 
Vertex(const CPoint &p) 
x = (Coord)p.x; y = (Coord)p.y; } 


operator CPoint () 

return CPoint((int)x, (int)y); } 

fendif 

Vertex& snap(Coord resolution) // Snap coords to given resolution 


if (resolution) 


{ 


x 
y. 


resolution * (int) (x/resolution) ; 
resolution * (int) (y/resolution) ; 


} 


return *this; 
}; 


class LineSegment; 
struct Rect 
{ 

Coord x0, yO, xl, yl; 


Rect (const Coord _x0 = 0, const Coord yO = 0, 
const Coord _xl = 0, const Coord yl = 0) 
: x0(_x0), yO(_y0), x1(_x1), yl(_yl) {} 


Rect (const Vertex &v0, const Vertex &vl): x0(v0.x), yO(vO.y), xl(vl1l.x), yl(vl.y) {} 
Rect (const Rect &r): x0(r.x0), yO(r.y0), xl(r.x1l), yl(r-yl) {} 


void normalize () 
{ 
if (xO > x1) 
swap (x0,x1); 
if (yO > yl) 
swap (y0,y1); 
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// Precondition: rectangle is normalized 
bool isVertexInside(const Vertex &v) const 
{ return v.x >= x0 && v.x <= xl && v.y >= yO && v.y <= yl; } 


// Precondition: rectangle is normalized 
void inflate(const Coord &dx, const Coord &dy) 
{ 


x0 -= dx; 
xl += dx; 
yO -= dy; 
yl += dy; 


}; 


class LineSegment 


{ 
static const float NOISE; // used for determining limit of parallel lines 


public: 

Vertex vertex0, vertexl; 

bool m_swappedVertices; // have vertices been swapped during sortVertices? 

Coord ymin, ymax; // y component of bounding box [don't use before call to 
sortVertices () ] 


typedef enum 
{ 


PARALLEL, // lines are parallel within tolerance level 
NO_INTERSECT, // lines segments don't intersect 
INTERSECT // line segments intersect 

} IntersectionType; 

LineSegment () {} 


// Sets bounding box of line segment to (vertex0.x, ymin, vertexl.x, ymax) 
void sortVertices () 
{ 
m_swappedVertices = !(vertex0 < vertexl); 
if (m_swappedVertices) 
swap (vertex0, vertexl1); 
if (vertex0.y < vertexl.y) 


ymin = vertex0.y; 
ymax = vertexl.y; 


else 


ymin = vertexl.y; 
ymax = vertex0.y; 


} 
LineSegment (const Vertex &vtx0, const Vertex &é&vtxl): vertex0(vtx0), vertexl(vtx1l) {} 


LineSegment (const LineSegment &1): vertex0(l.vertex0),vertexl(l.vertexl),ymin(l.ymin), 
ymax (l.ymax), m_swappedVertices(l.m_swappedVertices) {} 


IntersectionType calcIntersection(const LineSegment &l, Vertex &intersection, 
float é&alpha, float éalphal); 


bool isVertexNear(const Vertex &p, const Coord éresolution); 
// order by left edge of bounding box 


bool operator< (const LineSegment &ls) const 


{ 
return vertex0 < ls.vertex0; 


}; 


typedef list<Vertex>::iterator VtxIt; 
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as Intersection - to create the master list 

// Constructors 

// Intersection () 

// Intersection (const Intersection é&i ) 

// Intersection (const Vertex &vtx,LineSegment &li,LineSegment &1}j) 

// Intersection (const Vertex &vtx,Vertex é&o0rgl,Vertex s&org2,float &pl1, float 
&p2,LineSegment &li,LineSegment &1)) 


// Assignment 
Fy void operator= ( const Intersection éi) 


// Relational 
// bool operator< (const Intersection &i) const 
// bool operator<= (const Intersection &i) const 


al 


class Intersection 

{ 

public: 
Vertex v; 
float paraml,param2; 
int indexl, index2; 
Vertex originl, origin2; 
int winding; 
int direction; 
int selfIndex; 
LineSegment 11, 12; 


Intersection () {} 


Intersection (const Vertex &vtx, Vertex &o0rgl, Vertex ésorg2, 
float &pl,float &p2, LineSegment &li, LineSegment &1)) 


Vv = vtx; 
indexl = -1; 
index2 = -1; 


paraml = pl; 
param2 = p2; 
originl = orgl; 
origin2 = org2; 
11 = li; 
12 = 13; 

} 


Intersection ( const Vertex &vtx, LineSegment &li, LineSegment &1)j) 


originl = v; 

11 = li; 

12S 1g; 

winding = 0; 

direction = 0; 

selfIndex = 0; 
} 


Intersection( const Intersection &i): v(i.v), indexl(i.indexl), index2(i.index2), 
paraml (i.paraml),param2(i.param2), originl(i.originl), origin2(i.origin2), 
11(i.11), 12( 1.12 ), winding (i.winding), direction (i.direction), 
selfIndex(i.selfIndex) {} 


bool operator< (const Intersection &i) const 
{ return (v < i.v ); } 
bool operator<= (const Intersection &i) const 
{ return (v <= i.v ); } 


void operator= ( const Intersection &i) 
{ 
V=1.V; 
indexl = i.index1l; 
index2= i.index2; 
paraml = i.paraml; 
param2 = i.param2; 
originl = i.originl; 
origin2 = i.origin2; 
11 = i.11; 
12 = i.12; 
winding = i.winding; 
direction = i.direction; 


selfIndex = i.selfIndex; 


a2. 


typedef vector<Intersection> ::iterator InterIt; 


TITTILTTLTTLTITL TITTLE T TITAS 


// 
// 
// 
// 
// 
// 
// 
// 
// 
// 
// 
// 
// 


nVertex - temporary class to create the objects in the 
polygon edge array 

Constructors 

nVertex () 

nVertex ( const Vertex &vtx, float &p, int & i, LineSegment &12 ) 


nVertex ( const Vertex &vtx, int & i, LineSegment &12 ) 
nVertex ( const Vertex &vtx, LineSegment &12 ) 


Assignment 

void operator= ( const nVertex 6i) 
Relational 

bool operator< (const nVertex &n) const 


bool operator<= (const nVertex &i) const 


class nVertex 


{ 


public 


}; 


Vertex v; 
float param; 
int index; 
LineSegment 1; 
nVertex() {} 


nVertex ( const Vertex &vtx, float &p, int & i, LineSegment &12 ) 


nVertex ( const Vertex &vtx, int & i, LineSegment &12 ) 


vi = vtx; 
param = 0; 
index = i; 
1 = 12; 


nVertex ( const Vertex &vtx, LineSegment &12 ) 


vo = vtx; 
1 = 12; 
void operator= ( const nVertex 6i) 
v=i.v; 
1=i.1; 
param= i.param; 
index = i.index; 


} 


bool operator< (const nVertex &n) const 
{ return (param < n.param ); } 


bool operator<= (const nVertex &i) const 
{ return (param <= i.param ); } 


ae) 
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ae Pseudovertex - the polygon edge array 
// Constructors 

// Pseudovertex () 

// Pseudovertex (const nVertex &n) 

ST Pseudovertex (const Pseudovertex é&p) 

// 

// Relational 

Ty: bool operator< (const Pseudovertex &S) const 

// Comparison 

// bool operator== (const Pseudovertex &p) const 


class Pseudovertex 
{ 
public: 
vector<nVertex> ilist; 


Pseudovertex() {} 
Pseudovertex ( const nVertex én ) { ilist.push_back(n); } 
Pseudovertex ( const Pseudovertex &p ): ilist(p.ilist) {} 


bool operator< (const Pseudovertex &s) const 
{ return ( ilist.begin()->v < s.ilist.begin()->v ); } 


bool operator== (const Pseudovertex &p) const 
{ return ( ilist.begin()->v == p.ilist.begin()->v ); } 
}e 


typedef vector<Pseudovertex>::iterator PseudolIt; 


ZITTLLTTLTTLTITST TTI ATTA TITTIES 


// Poly: polygon; a cyclic collection of vertices 


// 

// Constructors 

// Poly() 

// Poly(const CPointé point): create triangle around point 

// Poly(const Poly& poly) 

// 

// Assignment 

// void operator= (const Polyé poly) 

// 

// Traversal 

ff Precondition: currPos must point to valid vertex. 

fa) Vertex& GetNextCirc(POSITION &currPos): for circular linkage 
It Vertex& GetPrevCirc (POSITION &currPos): for circular linkage 


class Poly 
{ 


public: 

list<Vertex> vtxList; 

Poly(const Poly &p): vtxList(p.vtxList) {} 

Poly() {} 

Poly(const Vertex &vtx, Coord resolution) 
vtxList.push_back(Vertex(vtx.x + 10.0f, vtx.y - 10.0f).snap(resolution) ); 
vtxList.push_back(Vertex(vtx.x - 10.0f, vtx.y - 10.0f).snap(resolution) ); 
vtxList.push_ back (Vertex (vtx.x , vtx.y + 10.0f).snap(resolution) ); 


void nextCirc(list<Vertex>::iterator &it) // for circular linkage 


if (++it == vtxList.end() 
ase = vtxList.begin(); 
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void prevCirc(list<Vertex>::iterator &it) // for circular linkage 


if (it == vtxList.begin () 
it = vtxList.end(); 
Sates 
} 
int size() const 


{ return vtxList.size(); } 


}; 
typedef list<Poly>::iterator PolyIt; 


TITITTTTTTTTTT TTT TTT TTT TTT AAA AAA AAA AAT AA AAA AAA AAA TTT TTT TTT TAA AAA AAA AAA AA AAA TTT 
// IntersectionList - list of polygon edges 
// Constructors 
// IntersectionList () 
// IntersectionList( const vector<Pseudovertex> é&p) 
class IntersectionList 
{ 
public: 
vector<Pseudovertex> p_list; 


IntersectionList () {} 


IntersectionList( const vector<Pseudovertex> &p) : p_list(p) {} 
}; 


TITITTTTTTTT TTT TTT TTT TATA AAA AAA AAA AAA AAA TA AAA AAT TTT TTT TAA AAA AAA AAA AAT 

// MultiPolygon: collection of polygons 

// 

TITTTTITTTTTTT TTT TTT TTT TATA AAA AAA AAA AAA AAA A AAA AAA ATTA TTT TTT TTA AAA AAA AAA AAT 

// Embedded class CMPPos: Position of vertex within multipolygon, (POSITION polygon, 
POSITION vertex) 


// 

// Constructors 

// CMPPos () 

// CMPPos(const POSITION _polyPos, const POSITION _vtxPos) 
// 

// Comparators 

// BOOL operator== (const CMPPos pos) const 

// BOOL operator!= (const CMPPos pos) const 

// 

// Assignment opertator 

[7 void operator= (const CMPPos pos) 

// 

If Predicate function to tell if position is valid 
// BOOL isValid() const 

// 


JITTILTTLITLTATLT TTT ATLL ATTA ATTA TITTTT A TTTTTTTT I ATTTTTT A TTATTTT 


// Constructors 


aie MultiPolygon() 

// MultiPolygon(const MultiPolygon& mPoly) 

// 

// Assignment 

Sf void operator= (const MultiPolygoné mPoly) 

// 

// Get position of multipolygon vertex within distance 'resolution' of a given vertex 
// MPPos GetActiveVtxPos (Vertex point, int resolution); 

// 


// For a given vertex. 'point', get position of multipolygon vertex vl 

// of a line (vl,v2) such that it is within a distance 'resolution' of the line. 
Sf MPPos GetActiveEdgePos (Vertex point, int resolution) 

// 
// Returns polygon containing vertex at given CMPPos position 
// Precondition: mpPos must be valid 

// Polygon& polygon(const CMPPos &mpPos) 

// 


// Returns vertex at given CMPPos position 


aye 


// Precondition: mpPos must be valid 


ais Vertex& vertex(const CMPPos &mpPos) 
class MultiPoly 
{ 
public: 
list<Poly> m_polyList; 
void random(Rect boundRect, int nVtx, int nPoly = 1); 
void dump (); 
bool read(const char* file name) ; 


bool write(const char* file name); 


bool writePS(const char* file name, bool bEOfill); 
struct MPPos 


{ 
list<Poly>::iterator m_polyIt; 
list<Vertex>::iterator m_vtxIt; 


MPPos(): m_vtxIt (NULL), m_polyIt(NULL) {} 
MPPos(PolyIt pi, VtxIt vi): m_polyIt(pi), m_vtxIt(vi) {} 
BOOL operator== (const MPPos pos) const 
return m_vtxIt == pos.m_vtxIt && m_polyIt == pos.m_polyIt; } 
BOOL operator!= (const MPPos pos) const 


return !(*this == pos); } 


void operator= (const MPPos pos) 


m vtxIt 
m_polyIt = 


= pos.m_vtxIt; 
pos.m polylIt; 


BOOL isValid() const 
return m_polyIt != 


NULL; } 
}; 


MultiPoly(const MultiPoly &mp): m_polyList(mp.m_polyList) 
MultiPoly() {} 
(const Vertex &vtx, 


MPPos GetActiveVtxPos 


MPPos GetActiveEdgePos(const Vertex &vtx, 


Rect boundingBox (void) ; 

vector<Vertex> findIntersections 
(IntersectionList &List, vector<Intersection> 

vector<Vertex> findMonotone (MultiPoly &resMPoly, 

void fillAddress(IntersectionList &ivList, vector 


void fillIndices(IntersectionList éivList, 


LineSegment::IntersectionType findIntersection( LineSegment &1l1, 


Vertex &intersection, float éalphal, float éalpha2); 


void gridify(Coord resolution); 


void polygonPartition(MultiPoly &resMPoly, 
vector<int> &windingVector) ; 


{} 


Coord resolution); 


Coord resolution); 


&tempList) ; 
vector<int> &windingVector) ; 


<Intersection> &interVector); 


vector <Intersection> &interVector) ; 


LineSegment &12, 


// Snap all vertex coordinates to given resolution 


vector<Intersection> &tempList, 


56 


// cross-product, pq x qr 


float xProd(const Vertex &p, Vertex &q, 


{ 
return ( (q.x * ( rey - p.y 
} 
}; 


typedef MultiPoly::MPPos MPPos; 
#endif 


)) 


ea 


(p.x* ( 


Vertex &r) 


q-y- r.y 


)) 


a7 


// Multipolygon.cpp 


#include "MultiPolygon.h" 

#include <cmath> 

#include <fstream> 

#include <algorithm> 

#include <iterator> 

#include <iostream> 

using namespace std; 

const float LineSegment::NOISE = le-5f; 


JITIISISTISISII AISI /// ine Intersection /////////////////1/1/1////1 
// 


// Calculates the point of intersection of two line segments. 

// 

// Precondition: both line segments have non-zero lengths. 

// Postcondition: The enumerated type function return value tells whether the lines 


// segments were parallel, non-intersecting, or intersecting. 
// "intersectionPoint" is valid iff the returned type is not PARALLEL. 
// 


LineSegment::IntersectionType LineSegment::calcIntersection 
(const LineSegment &1, Vertex &intersection, float &alpha, float éalphal) 


{ 


float dx21 = vertexl.x - vertex0.x, dy21 = vertexl.y - vertex0.y, 
dx43 = l.vertexl.x - l.vertex0.x, dy43 = l.vertexl.y - l.vertex0.y; 
float dem = dx21 * dy43 - dy21 * dx43; 


// parallel lines 

if (fabs(dem) < NOISE) 
return PARALLEL; 

else { 
float dx1l2 = vertex0.x - l.vertex0.x, 

dyl2 = vertex0.y - l.vertex0.y; 

float dx = l.vertex0.x - vertex0.x, 

dy = l.vertex0.y - vertex0.y; 


alpha = (dyl2*dx43 - dx12*dy43) /dem; 


alphal = ( dy21*dx - dx21*dy) /dem; 
// The intersecting point 


intersection.x = vertex0.x + dx21 * alpha; 
intersection.y = vertex0.y + dy21 * alpha; 


// test for segment intersecting (alpha) 


if ((alpha < 0.0) || (alpha > 1.0) 
return NO_ INTERSECT; 
else { 
float num = dyl2*dx21 - dx12*dy21; 
if (dem > 0.0) { 
if (num < 0.0 || num > dem) 
return NO INTERSECT; 


if (num > 0.0 || num < dem) 
return NO INTERSECT; 


} 
} 
return INTERSECT; 


} 


bool LineSegment::isVertexNear (const Vertex &vtx, const Coord &resolution) 

/* 

* Determine whether point vtx is within distance 'resolution' from line segment. 
* 
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* Written by T. Hain, Jan 1998. 

Af 

{ 
// CRect activeRect (vertex0.m_point, vertexl.m point); 
// activeRect.NormalizeRect (); 


// activeRect.InflateRect (resolution, resolution); // for vertical or horizontal 
edges 

// if (l!activeRect.PtInRect (*this) ) 

iy return false; // trivial reject (i.e., *this is not within line seg's 


bounding box) 
//Rectangle rect (*this) ; 
//cect.inflate (resolution, resolution) ; 
//if (!rect.isVertexInside (vtx) ) 
// veturn false; 


Coord xl = vertex0.x, 
yl = vertex0.y; 
Coord x2 = vertexl.x, 
y2 = vertexl.y; 
Coord a = (yl - y2), 
by SCS &2)5 
ee Steal & ye soy) sey # x2) a1) 


Coord xba = x2 - xl, 

yba = y2 - yl; 

double lsqr = xba * xba + yba * yba, 
1 = sqrt(lsqr); 

Coord yac = yl - vtx.y, 


xac = xl - vtx.x; 
double r = ( -yac * yba - xac * xba); 
double s = ( yac * xba - xac * yba) / 1; 


return (fabs(s) < resolution && r > -resolution && r < lsqr + resolution); 


} 


TITISIITITIIIIII IIIT MaltiPoly J//IIIIIIIIIIIIII11/// 
// 
MultiPoly::MPPos MultiPoly: :GetActiveVtxPos (const Vertex &vtx, Coord resolution) 
// Find first vertex such that vtx is within a distance 'resolution' of it. 
// Output MPPos for that 'active' vertex, or invalid if none is found 
{ 

Rect activeArea(vtx, vtx); 

activeArea.inflate(resolution, resolution); 

for (PolyIt pi = m polyList.begin(); pi != m_polyList.end(); ++pi) 

for (VtxIt vi = pi->vtxList.begin(); vi != pi->vtxList.end(); ++vi) 
if (activeArea.isVertexInside (*vi) ) 
return MPPos(pi,vi); 

return MPPos (NULL, NULL) ; 

} 


MultiPoly::MPPos MultiPoly: :GetActiveEdgePos (const Vertex &vtx, Coord resolution) 
// Find first edge such that vtx is within a distance 'resolution' of it. 
// Output (activePolyIter,activeVtxIter) is only valid if active edge is found, and 
// vefers to vertex at one end of edge. Vertex at other end of edge is 
// (activePolyPos, GetNextCirc (activeVtxPos) ). 
// Returns true if an active edge is found 
{ 
for (PolyIt pi = m polyList.begin(); pi != m_polyList.end(); ++pi) 
{ 
VtxIt vi = pi->vtxList.begin(); 
Vertex vO = *vi; 
for (int i = 0; i < pi->size(); ++i) 
{ 
pi->prevCirc(vi); 
Vertex vl = *vi; 
LineSegment 1(v0,v1); 
if (l.isVertexNear(vtx, resolution) ) 
return MPPos(pi,vi); 
vO = vil; 
} 
} 
return MPPos (NULL, NULL) ; 
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} 


bool MultiPoly::writePS(const char *file name, bool bEOfill) 
{ 

ofstream out (file name); 

Rect rect (boundingBox()); 

out << "/poly {newpath"; 


for (PolyIt pi = m polyList.begin(); pi != m_polyList.end(); ++pi) 
{ 
out << "\n\t" << pi->vtxList.begin()->x << " " << pi->vtxList.begin()->y 
<< " moveto\n"; 
for (VtxIt vi = pi->vtxList.begin()++; vi != pi->vtxList.end(); ++vi) 
out << "\t" << vi->x << " " << vi->y << " lineto\n"; 


out << "\tclosepath\n"; 

} 

out << "} def\n\t16 774 translate\n\t0.8 setgray\n\tpoly " 
<< (bEOfI11>? “eofill\n" «+ “ELIIA\n") 
<< "\t0O setgray\n\tpoly stroke\n\tshowpage"; 


return true; 


} 


bool MultiPoly::write(const char *file name) 
{ 

ofstream out (file name); 

out << m polyList.size() << "\n"; 


for (PolyIt pi = m_polyList.begin(); pi != m_polyList.end(); ++pi) 
{ 
out << pi->vtxList.size() << ""; 
for (VtxIt vi = pi->vtxList.begin(); vi != pi->vtxList.end(); ++vi) 
out << vi->y << " " << vi->x << ""; 


} 
out: 220 as 
return true; 


} 


bool MultiPoly::read(const char *file name) 
{ 
ifstream in(file name); 
m polyList.clear(); 
list<Poly>::size type nPoly(*istream_iterator<typedef list<Poly>::size type>(in)); 
for (list<Poly>::size type p = 0; p < nPoly; ++p) 
{ 
Poly poly; 
list<Vertex>::size type nVtx(*istream_iterator<list<Vertex>::size type>(in)); 
for (list<Vertex>::size type v = 0; v < nVtx; ++v) 
{ 
Vertex vtx(*istream_iterator<Coord>(in),*istream_iterator<Coord>(in)); 
poly.vtxList.push_back(vtx); 
} 
m_polyList.push_ back (poly); 
} 
return true; 


} 


void MultiPoly::dump () 
{ 
ofstream out ("dump.txt"); 
out << m_polyList.size() << "\n"; 


for (PolyIt pi = m_polyList.begin(); pi != m_polyList.end(); ++pi) 
{ 
out << pi->vtxList.size() << ""; 
for (VtxIt vi = pi->vtxList.begin(); vi != pi->vtxList.end(); ++wvi) 
out << vi->y << " " << vi->x << ""; 


} 


out << "\n"; 
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// generates a random multi polygon within the given rectangular area, 
// for a given number of vertices and number of polygons 
void MultiPoly::random(Rect boundRect, int nVtx, int nPoly) 
{ 
srand( (unsigned)time( NULL ) ); 
boundRect.normalize(); 
m_ polyList.clear(); 
for (int p = 0; p < nPoly; ++p) 
{ 
Poly poly; 
for (int v = 0; v < nVtx; ++v) 
{ 
Vertex vtx 
( 
boundRect.x0 + rand()%((int) (boundRect.xl-boundRect.x0)), 
boundRect.y0 + rand()%((int) (boundRect.yl-boundRect.y0) ) 
)e 
poly.vtxList.push_back(vtx) ; 
} 
m_polyList.push_ back (poly); 


} 


Rect MultiPoly: :boundingBox () 
{ 
Rect rect (*m_polyList.begin()->vtxList.begin(), *m_polyList.begin()->vtxList.begin()); 
for (PolyIt pi = m polyList.begin(); pi != m_polyList.end(); ++pi) 
for (VtxIt vi = pi->vtxList.begin(); vi != pi->vtxList.end(); ++vi) 
{ 


rect.x0, vi->x); 


rect.x0 = min ) 

rect.y0, vi->y); 
) 
) 


rect.y0 = min 
rect.xl = max 
rect.yl = max 


, 


rect.xl, vi->x 
rect.yl, vi->y 


a 


} 


return rect; 


} 


// finds the intersection point between two line segments and thier parametric 
// values with respect to each intersecting line segment 
LineSegment::IntersectionType MultiPoly::findIntersection( LineSegment &11, 
LineSegment &12, Vertex &intersection, float &alphal, float &alpha2) 

{ 
LineSegment 13(11); 

LineSegment 14(12); 

1l.sortVertices (); 

12.sortVertices (); 

Vertex vtx; 

LineSegment::IntersectionType resultType; 

if ( 1l.vertexl < 12.vertex0 ) // no x-overlap 


// check for y-overlap 
te £2) 12) 
{ 
if (ll.ymax < 12.ymin ) 
return (LineSegment::NO_ INTERSECT ); 


if (12.ymin < 11.ymax) 
return LineSegment::NO_INTERSECT; 
} 
else 
// check for adjacency 
if ( (1ll.vertex0 == 12.vertex0) || (11l.vertexl == 12.vertexl) || 
(ll.vertexO == 12.vertexl) || (1ll.vertexl == 12.vertex0) ) 
return ( LineSegment::NO_ INTERSECT ); 
resultType = 13.calcIntersection(14, vtx, alphal, alpha2 ); 
if ( resultType == LineSegment::INTERSECT ) 
{ 


intersection = vtx; 
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return resultType; 
} 
return (resultType); 


} 


// this is the method called to fill the polygon useing NZW algorithm 
// finds the intersections using trivial rejections considering the 
// bounding box overlap and monotonic chains 


vector<Vertex> MultiPoly::findMonotone (MultiPoly &resMPoly, vector<int> &windingVector) 


{ 
vector<Vertex> intersections; 
vector<Intersection> tempList; 
vector<Pseudovertex> p_list; 
p_list.clear(); 
// add the vertices to the ivList in the form of Intersection 


for ( PolyIt pii = m_polyList.begin(); pii != m_polyList.end(); ++pii) 


{ 
VtxIt vit0O, vitl; 


for (VtxIt vi = pii->vtxList.begin(); vi != pii->vtxList.end();++vi) 


{ 

Pseudovertex pseudo; 

vitO = vi; 

vitl = vi; 
pii->prevCirc (vit0); 
pii->nextCirc(vitl); 
LineSegment 11(*vit0,*vi); 
LineSegment 12(*vi, *vitl); 
12.sortVertices (); 
nVertex i(*vi, 12); 
tempList.push_ back(Intersection ( *vi, 12, 11) ); 
pseudo.ilist.push_ back (i); 
p_list.push_back( pseudo ); 


} 
IntersectionList ivList(p_ list); 


// add the intersection points to the list, as Intersection objects 
for (PolyIt pi = m polyList.begin(); pi != m_polyList.end(); ++pi) 
{ 

VtxIt vi0 = pi->vtxList.begin(), vil; 

unsigned int x _monotone=0; 

bool reached = false; 


vil = vi0; 

VtxIt start = vi; 

++vil; 

//find intersection points for every polygon. 

while ( vil != pi->vtxList.end() && (reached == false) ) 


{ 


LineSegment 11(*vi0, *vil); 


vid = vil; 
++vil; 
if ( vil == pi->vtxList.end() ) 


vil = start; 
reached = true ; 
} 
LineSegment 12(*vi0, *vil); 

LineSegment temp(12 ); 

VtxIt walker, walkerl; 

walker = start; 

walkerl = walker; 

++walkerl; 

// walk thru the list from the beginning to find 
// intersection with the non-monotone edge 

while ( walkerl != vi0 ) 


LineSegment 11(*walker, *walkerl ); 
12 = temp; 
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Vertex intersection; 

float all, al2; 

if ( findIntersection(11,12,intersection, all, al2) 
== LineSegment::INTERSECT ) 


intersections.push back(intersection) ; 
Vertex searchVertexl, searchVertex2; 
if (11.m_swappedVertices) 


searchVertexl = 11.vertexl; 
else 
searchVertexl = 11.vertex0; 
if ( 12.m_swappedVertices) 
searchVertex2 = 12.vertexl; 
else 
searchVertex2 = 12.vertex0; 
tempList.push_back( Intersection ( intersection, searchVertex1l, 
searchVertex2, all, al2, 11, 12 ) ); 


} 
walker= walkerl; 
++walkerl; 


} 
vector<Vertex> resultIntersections = findIntersections(ivList, tempList); 


for ( unsigned int i=0; i < resultIntersections.size() ; it+ ) 
intersections.push_ back(resultIntersections[i]); 


sort ( ivList.p list.begin(), ivList.p list.end() ); 


fillAddress(ivList, tempList); 
polygonPartition( resMPoly,tempList,windingVector) ; 
return intersections; 


} 


// finds the intersection points between every polygon edge 
vector<Vertex> MultiPoly::findIntersections(IntersectionList é&ivList, 
vector<Intersection> &tempList) 
{ 
vector<LineSegment> edges, edges1; 
vector<Vertex> intersections; 
for (PolyIt pi = m polyList.begin(); pi != m_polyList.end(); ++pi) 
{ 
PolyIt pid = pi; 
VtxIt vi0 = pi0->vtxList.begin(), vil; 
unsigned nVtx = pi0->vtxList.size(); 
for (unsigned i = 0; i < nVtx; ++i) 
{ 
vil = vi; 
pi0->nextCirc (vil); 
Vertex v0(*vi0),vl1(*vil); 
LineSegment edge(*vi0, *vil); 
edge.sortVertices(); 
edges.push_back (edge) ; 
vid = vil; 
} 
sort (edges.begin(), edges.end()); 
++pi0; 
while ( piO != m_ polyList.end() ) 
{ 
VtxIt vi2 = pi0->vtxList.begin(), vi3; 
unsigned nVtx = pi0->vtxList.size(); 
for (unsigned i = 0; i < nVtx; ++i) 
{ 
vi3 = vi2; 
pi0->nextCirc(vi3); 
Vertex v0 (*vi2),vl1(*vi3); 
LineSegment edge(*vi2, *vi3); 
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edge.sortVertices(); 
edges1l.push_back (edge) ; 
vi2 = vi3; 
} 
sort (edgesl.begin(), edgesl.end()); 
// make vector of intersections 
unsigned nEdges = edges.size(); 
unsigned nEdgesl = edgesl.size(); 
if (nEdges > 2) 
{ 
for (unsigned i = 0; i < nEdges; ++i) 
for (unsigned j = 0; j < nEdgesl; ++3) 
{ 
LineSegment li, lj; 


if ( !edges[i].m_swappedVertices ) 
li = edges[i]; 

else 
li.vertex0 = edges[i].vertexl1; 
li.vertexl = edges[i].vertex0; 

if ( !edgesl[j].m_swappedVertices ) 
1j = edgesl[j]; 

else 
1j.vertex0 = edgesl[j].vertexl; 
1j.vertexl = edgesl[j].vertex0; 


// Because of sort above: edges[i].rect.x0 <= edges[j].rect.x0 


// if bounding boxes have x-overlap 
if (edgesl[j].vertex0.x < edges[i].vertexl.x) 
{ 
// check for y-overlap of bounding boxes 
if (edges[i].ymin < edgesl[j].ymin) 
{ 
if (edges[i].ymax <= edgesl[j].ymin) continue; 
} 
else 
if (edgesl[j].ymax <= edges[i].ymin) continue; 
// Check if both left [right] ends are not coincident (i.e., 
// valid intersection) 
if (edges[i].vertex0 != edgesl[j].vertex0 && edges[i].vertexl 
!'= edges1[j].vertex1) 


Vertex intersection; 

float paraml, param2; 

if (li.calcIntersection(lj, intersection, paraml, param2) 
== LineSegment: : INTERSECT) 


intersections.push_ back(intersection) ; 


Vertex searchVertexl, searchVertex2; 
if (edges[i].m_swappedVertices) 
searchVertexl = edges[i].vertexl; 
else 
searchVertexl = edges[i].vertex0; 
if ( edgesl[j].m_swappedVertices) 
searchVertex2 = edges1l[j].vertex1; 
else 
searchVertex2 = edgesl[j].vertex0; 
li.sortVertices (); 
14j.sortVertices(); 
tempList.push_back( Intersection ( intersection, 
searchVertexl, searchVertex2, paraml, param2, 


is ea Og) 
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} 
++pi0; 
}// end of while for pi0 counter 
}// end of for 
return intersections; 


} 


// a function to fill polygon edge array(ivList) with the indices from the 
// intersection master list 
void MultiPoly::fillAddress(IntersectionList &ivList, vector<Intersection> é&interVector) 
{ 
sort ( interVector.begin() , interVector.end() ); 
for ( int i =0; i < interVector.size() ; i ++ ) 
{ 
if ( (interVector[i].paraml != 0 ) && (interVector[i].param2 != 0 ) ) 
{ 
int ill = 0, i12 = 0, i21 = 0, i22 = 0; 
Vertex searchVertexl = interVector[i].originl; 
Vertex searchVertex2 = interVector[i].origin2; 
WV 


while ( ivList.p_ list[ill].ilist[0]. !'= searchVertex1) 
illt++; 

while ( ivList.p_ list[i21].ilist[0].v != searchVertex2) 
i21++; 


nVertex Interl(interVector[i].v, interVector[i].paraml, i, 
interVector[i].11), Inter2(interVector[i].v, interVector[i].param2, i , 
interVector[i].12); 

ivList.p list[ill].ilist.push_back(Inter1); 

ivList.p list[i21].ilist.push_back(Inter2) ; 

} 
else 
{ 
int vtxIt = 0; 


Vertex vtx = interVector[i].v; 
while ( ivList.p_list[vtxIt].ilist[0].v != vtx ) 
vtxItt+t+; 


ivList.p list[vtxIt].ilist[0].index = i; 


} 

for ( PseudoIt ps = ivList.p_list.begin(); ps != ivList.p list.end(); ++ps) 
sort ( ps->ilist.begin() + 1, ps->ilist.end() Yaz 

fillIndices( ivList, interVector ); 


} 


// A method to fill in the last node of polygon edge array (ivList), which points to 
// the ending vertex of the line 
void MultiPoly::fillIndices ( IntersectionList &ivList, vector<Intersection> 
é&interVector) 
{ 
// setting up the last Vertex that points to the starting vertex of the 
// intersecting and setting up its corresponding index in the intersection 
// master list (interVector) 
for ( PseudoIt ps = ivList.p_list.begin(); ps != ivList.p list.end(); ++ps) 
{ 
int vtxIt = 0; 
Vertex searchVtx; 
LineSegment ls( (ps->ilist.end() - 1)->1 ); 
if ( 1ls.m_swappedVertices) 
searchVtx = 1ls.vertex0; 
else 
searchVtx = ls.vertexl; 


while ( interVector[vtxIt].v != searchVtx ) 
Vex Gah; 
ps->ilist.push_back( nVertex( searchVtx, vtxIt, (ps->ilist.end() - 1)->1l yy ) og 
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// setting up the indices for the intersections in the interVector 
for ( ps = ivList.p list.begin(); ps != ivList.p list.end(); ++ps) 


for (int i= 0; i < ps->ilist.size()-1l; ++i) 


{ 
nVertex tempVtx = ps->ilist[itl]; 
int tempIndex = ps->ilist[i].index; 


if ( interVector[tempIndex].origin2 == ps->ilist[0].v ) 
interVector[tempIndex].index2 = tempVtx.index; 


else 
if ( interVector[tempIndex].originl == ps->ilist[0].v ) 
interVector[tempIndex].indexl = tempVtx. index; 
else 
if ( interVector[tempIndex].v == ps->ilist[0].v ) 


interVector[tempIndex].indexl = tempVtx.index; 


} 

// setting up the self-index of each interseciton object 

for ( int i =0 ; i < interVector.size() ; i ++ ) 
interVector[i].selfIndex = i; 


// Snap all vertex coordinates to given resolution 
void MultiPoly::gridify(Coord resolution) 
{ 


if (resolution != 0) 
for (PolyIt pi = m polyList.begin(); pi != m_polyList.end(); ++pi) 
for (VtxIt vi = pi->vtxList.begin(); vi != pi->vtxList.end(); ++vi) 


vi->snap (resolution) ; 


// traverse from the left most vertex of the ivList and form the queue of subsequent 
// polygons done until the queue becomes empty i.e all the vertices of the 
// multipolygon has been processed 
void MultiPoly::polygonPartition(MultiPoly &resMPoly, vector<Intersection> &tempList, 
vector<int> é&windingVector) 
{ 
// vesulting polygon 
vector<Intersection> interVector; 
Polylt pile; 
VtxIt vit; 
bool finished = false; 


Vertex startVtx, currVtx; 

Intersection previnter, temp; 

interVector = tempList; 

int i =0, currWinding = 0, startIndex = 0, firstIndex = 0, secondIndex = 0; 

bool foundWinding = false; 

int currEO =0; 

stack<Intersection> interStack; 

MultiPoly resMP; 

vector<int> wVector; 

int currDirection =0; 

do 

{ 
int index =0, currIndex;// index of interVector 
Poly currPoly; 


if ( !interStack.empty() ) 

{ 
while ( !interStack.empty() ) 
{ 


index = interStack.top().selfIndex; 
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temp = interStack.top(); 
if ( interVector[index].indexl == -1 && interVector[index].index2 == -1) 
{ 
interStack.pop(); 
} 


else 
break; 
} 
} 
if ( !interStack.empty() && interStack.top().selfIndex > 0) 
index= interStack.top().selfIndex; 
else 
{ 
index=0; 
while(( interVector[index].indexl == -1) 
&& ( interVector[index].index2 == -1 ) &&(index < interVector.size()) ) 


{ 
temp = interVector [index]; 
index++; 


} 

if ( index >= interVector.size()-1 
finished = true; 

else 
currIndex = index; 


if ( !finished ) 
{ 


startVtx = interVector[currIndex] .v; 


currVtx = interVector[currIndex].v; 
previnter = interVector[currIndex]; 
startIndex = currIndex; 


firstIndex = startIndex; 
Vertex vtx; 


// fix vertex to find the winding number before changing the index values!! 


if ( interVector[startIndex].indexl == -1 ) 
{ 
if ( interVector[startIndex].11.m_swappedVertices) 
vtx = interVector[startIndex].1ll.vertex1; 
else 
vtx = interVector[startIndex] .ll.vertex0; 


if ( interVector[startIndex] .12.m_swappedVertices) 
vtx = interVector[startIndex] .12.vertex1; 

else 
vtx = interVector[startIndex] .12.vertex0; 


} 
if ( interVector[startIndex].indexl == -1 || interVector[startIndex] .index2 
ae et 


currWinding = interVector[startIndex] .winding; 


currPoly.vtxList.push_back( currVtx); 
interStack.push (interVector[currIndex]); 


do { 
if (( interVector[currIndex].indexl == -1) 
&& ( interVector[currIndex].index2 == -1 ) ) 
currIndext+; 
else 
{ 
if ( ( interVector[currIndex].originl == currVtx ) 
&& ( interVector[currIndex].index2 == -1 ) ) 


int tempIndex = interVector[currIndex].index1; 
previnter = interVector[currIndex]; 
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interVector[currIndex].indexl = -1; 
currIndex = tempIndex; 
} 
else 
{ 
if ( (( interVector[currIndex].originl == prevInter.originl) 
|| ( interVector[currIndex].originl == prevInter.origin2 )) 
&& ( interVector[currIndex].index2 != -1 ) ) 
{ 
int tempIndex = interVector[currIndex] .index2; 
previnter = interVector[currIndex]; 
interVector[currIndex].index2 = -1; 
currIndex = tempIndex; 
} 
else 


{ 
int tempIndex = interVector[currIndex].index1; 
previnter = interVector[currIndex]; 
interVector[currIndex].indexl = -1; 
currIndex = tempIndex; 


(currIndex > -1) 


currVtx = interVector[currIndex] .v; 
currPoly.vtxList.push_back( currVtx); 
interStack.push (interVector[currIndex]) ; 


if ( !foundWinding) 
{ 
if (firstIndex == startIndex) 
firstIndex = currIndex; 
else 
{ 
secondIndex = currIndex; 
foundWinding = true; 


Intersection templ,temp2; 


templ = interVector[firstIndex]; 
temp2 = interVector[secondIndex]; 
if (( interVector[startIndex].v.x 
== interVector[startIndex].originl.x) 
&& ( interVector[startIndex].v.y 
== interVector[startIndex].originl.y ) 
|| ( interVector[startIndex] .v.x 
== interVector[startIndex].origin2.x) 
&& ( interVector[startIndex].v.y 
== interVector[startIndex].origin2.y )) 
{ 
float winl; 
winl = xProd(interVector[startIndex] .12.vertex0, 
interVector[startIndex].v, interVector[firstIndex].v); 
if (winl > 0 
currDirection = 0; 
currWinding--; 
else 
currDirection = 1; 
currWindingt+; 


interVector[startIndex].direction = currDirection; 
} 
else 
{ 
float win; 
win = xProd(vtx, interVector[startIndex].v, 
interVector[firstIndex].v); 


if (win > 0 ) 
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} 


bo 
while 


} 


wvector. 


currDirection = 0; 


if ( interVector[startIndex].direction != currDirection) 
currWinding= currWinding - 2; 
else 
currWinding--; 
} 
else 
{ 
currDirection = 1; 
if ( interVector[startIndex].direction != currDirection) 
currWinding= currWinding + 2; 
else 
currWindingt+; 
} 
} 
interVector[startIndex].direction = currDirection; 
interVector[firstIndex].direction = currDirection; 
interVector[secondIndex].direction = currDirection; 
interVector[startIndex] .winding = currWinding; 
interVector[firstIndex] .winding = currWinding; 
interVector[secondIndex] .winding = currWinding; 


else 


} 


while ( (currVtx 


interVector[currIndex].winding = currWinding; 
interVector[currIndex].direction = currDirection; 


'= startVtx) && ( currIndex > -1 ) ) ; 


push back (currWinding) ; 


foundWinding = false; 


resMP.m_ 


( 


resMPoly 
windingVector = wVector; 


!'finished ); 
= resMP; 


polyList.push_back(currPoly) ; 
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// SLA.h 
#include "MultiPolygon.h" 
#include <fstream> 


using namespace std; 

TITITTTTTTT TTT TTT TTT TATA ATTA AAA AAA AAA AAA AAA AAA AAA ATT TTT TTT TATA AAA AAA AA AAA AAT 
// Global Edge Table entry 

// Constructors 

// GET() 

// GET ( const LineSegment &1, Coord yl, Coord yO, Coord x0, Coord m ) 

// GET( const GET &g): 11(g.11), ymin(g.ymin), ymax (g.ymax), 


// xOfymin(g.xOfymin) , slope(g.slope) 
// Assignment 
// void operator= ( const GET &g) 


// Relational 
// wool operator< (const GET &g) const 
// bool operator<= (const GET &g) const 


class GET 

{ 

public: 
LineSegment 11; 


public: 
Coord ymax, ymin, xOfymin, slope; 


GET() {} 


GET ( const LineSegment &1, Coord yl, Coord y0, Coord x0, Coord m ) 
{ 


11 = 1; 
ymin = y0; 
ymax = yl; 


slope = m; 
xOfymin = x0; 


} 


GET( const GET &g): 1l1(g.11), ymin(g.ymin), ymax(g.ymax), 
xOfymin(g.xOfymin),slope(g.slope) {} 


bool operator< (const GET &g) const 
{ return (ymin < g.ymin ); } 


bool operator<= (const GET &g) const 
{ return (ymin <= g.ymin ); } 


void operator= ( const GET &g) 
{ 

11 = g.11; 

ymin = ymin; 
ymax = ymax; 
xOfymin = g.xOfymin; 
slope = g.slope; 


1 
g 
g 


}; 


TITILITTTLSTTTTALS TTT T TTT TTT TTT TTA TT TTA TTT TAT ATT TATA TTT TAA TTT TT 
// Active Edge Table Entry 

// Constructors 

// AETEntry () 

// AETEntry (const LineSegment &l, float y, float x, float m_inv ) 
Ms AETEntry (const AETEntry &a) 

// Assignment 

// wool operator= ( AETEntry *aet) 

// Relational 

// bool operator<= (const AETEntry &aet) const 

// bool operator< (const AETEntry é&aet) const 
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class AETEntry 
{ 
public: LineSegment edge; 
float ymax, xval,m_inverse; 


AETEntry () {} 


AETEntry ( const LineSegment &l, float y, float x, float m_inv ) 
{ 


edge = 1; 
ymax = y; 
xval = x; 
m_inverse = m_inv; 


} 


AETEntry ( const AETEntry &a ) 
edge (a.edge ), ymax(a.ymax), xval(a.xval), m_inverse(a.m_inverse) {} 
bool operator<= (const AETEntry é&aet) const 
{ 
return (xval <= aet.xval ) ; 


} 


bool operator= ( AETEntry *aet) 
{ 

edge = aet->edge; 

ymax = aet->ymax; 

xval = aet->xval; 

m_inverse = aet->m_inverse; 


} 


bool operator< (const AETEntry &aet) const 
{ return (xval < aet.xval ); } 
}; 


typedef vector<LineSegment> edges; 

typedef vector<GET> GETable; 

typedef vector<GET>::iterator get_it; 
typedef vector<AETEntry> AET; 

typedef vector<AETEntry>::iterator aet_it; 


TITILITTTTITTTTATS TATA TTT TTT ATTA ATTA TTT ATTA TTT TTA TTT ATTA 

// SLAM - creates an SLA objetc for a multi polygon 

// contains the Global Edge Table and Active Edge Table 

// Constructors 

// SLAM () 

// SLAM ( MultiPoly &m ) 

// SLAM ( const GETable &get, AET &aet ) 

class SLAM 

{ 

public: 
MultiPoly mp; 
GETable polyGET; 
BET polyAET; 


SLAM () {} 


SLAM ( MultiPoly é&m ) 
{ mp = m; } 


SLAM ( const GETable &get, AET &aet ) 
{ polyGET = get; polyAET = aet; } 


void DrawSpan (float &y, AETEntry *pl, AETEntry *p2, CDC &pDC, bool isEO, 
int &winding, bool even) ; 

void FillPoly (MultiPoly *mpoly, CDC &pDC, bool isEO ); 

vector <LineSegment> GetEdges( MultiPoly *m_poly , float &ymax) ; 

void InitGET(GETable &get, vector<LineSegment> &edgesVector , float y); 

void UpdateAET(AET &aet, float y, GETable &get); 

void UpdateAETEntry(AETEntry &g); 
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void UpdateGET( GETable é&get, float &y); 


float xProd(const Vertex &p, Vertex &q, Vertex &r) 
(return: <(iqesee Cory ape yi eqs Ss ary op aM RA OC pay sigh yi ai oF 
}; 
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// SLA.cpp 
#include "SLA.h" 
#include <iostream> 
using namespace std; 
/* draws a span between two Active Edge Table entries of the polygon */ 
void SLAM::DrawSpan (float &y, AETEntry *pl, AETEntry *p2, CDC &pDC, 
bool isEO,int &éwinding, bool even) 

{ 

float xl, x2; 

/* don't draw spans with exactly one integer point */ 

if (pl == p2 

return; 
if ( pl->xval > p2->xval ) 


xl = p2->xval; 
x2 = pl->xval; 


else 


xl = pl->xval; 
x2 = p2->xval; 


/* draw pixels */ 
if (isEO && even) 


if (( xl >0) && ( x2 > 0) && (xl <= max(pl->xval,p2->xval) ) 
&& ( x2 <= max(pl->xval,p2->xval ) )) 
for (float x = xl; x <= x2; x+t+) 
pDC.SetPixelV( (int)x, (int) y, RGB(213, 175,255 )); 


if ( !isEO ) 


// find the direction of the segment 
if ( pl->edge.vertex0 > pl->edge.vertexl ) 


windingt+t+; 
else 
winding--; 
if ( winding != 0 ) 
if (( xl >0) && ( x2 > 0) && (xl <= max(pl->xval,p2->xval) ) 


&& ( x2 <= max(pl->xval,p2->xval ) )) 
for (float x = xl; x <= x2; x+t+) 
pbC.SetPixelV((int) x, (int)y, RGB(213, 175,255 )); 
if ( p2->edge.vertex0 > p2->edge.vertexl ) 
windingt+t+; 
else 
winding--; 


} 


// fill an arbitrary polygon with the even-odd interior rule 
void SLAM :: FillPoly (MultiPoly *mpoly, CDC &pDC, bool isEO) 
{ 

float y = 0; 

bool even = true; 

GETable get; 

AET aet; 

float ymax = 0; 

int winding ; 

edges eVector; 

// ScanBuffer myScanBuff; 

aet_it pl,p2,p; 

eVector = GetEdges(mpoly, ymax); 

//initialise scanbuffer; 

/* create Global Edge Table, initialize y and AET */ 

InitGET (get, eVector,y); 

y =get.begin()->ymin; 

/* sweep the scanline */ 

while (y < ymax) 
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/* insert / delete edges in Active Edge Table */ 
UpdateAET (aet, y, get); 
// insert or delete edges from get 
if (!get.empty () 
UpdateGET (get,y); 
/* sort Active Edge Table by increasing x */ 
sort (aet.begin(), aet.end()); 
if ( !aet.empty() ) 
if ( aet.begin()->edge.vertex0 > aet.begin()->edge.vertexl ) 
winding =1; 
else 
winding = -1; 
/* draw spans with even-odd rule */ 
even = true; 


for ( pl = aet.begin() ;pl != aet.end(); pltt) { 
p2 = pl; 
p2tt+; 
if ( p2 != aet.end() ) 


{ 
// DrawSpan (y, pl, p2, pDC,isEO, winding,even ); 
DrawSpan (y, (AETEntry *)&(*pl), 
(AETEntry *)&(*p2), pDC,isEO,winding,even ); 
} 
even = !even; 
} 
/* compute intersections for next scanline */ 
ytt; 
for (p = aet.begin(); p != aet.end(); ptt ) 
UpdateAETEntry (*p); 
} 
} 
// Convert the given multi polygon into a vector of polygon edges 
// and also keeps track of the ymax of the given polygon 
vector<LineSegment> SLAM::GetEdges( MultiPoly *m_poly, float &ymax ) 
{ 
vector<LineSegment> edges1l; 
PolyIt pi = m_poly->m_polyList.begin(); 
// for every polygon 
PolyIt pid = pi; 
// get the polygon vertex list 
VtxIt vid = pi0->vtxList.begin(), vil; 
unsigned nVtx = pi0->vtxList.size(); 
for (unsigned i = 0; i < nVtx; ++i) 
{ 
vil = vi0; 
pi0->nextCirc (vil); 
Vertex v0(*vi0),vl1l(*vil); 
// set ymax 
if ( ( vi0->y > ymax )&&(vi0->y > vil->y ) ) 
ymax = vi0->y; 
else 
if ( vil->y > ymax ) 
ymax = vil->y; 
LineSegment edge(*vi0, *vil); 
edgesl.push_back (edge) ; 
vid = vil; 
} 
// get the next polygon 
++pi0; 
while ( piO != m_poly->m_polyList.end() ) 
{ 
VtxIt vi2 = pi0->vtxList.begin(), vi3; 
unsigned nVtx = pi0->vtxList.size(); 
for (unsigned i = 0; i < nVtx; ++i) 
{ 
vi3 = vi2; 
pi0->nextCirc(vi3); 
Vertex v0 (*vi2),vl1(*vi3); 
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if ( (wi2->y > vi3->y) && ( vi2->y > ymax ) ) 
ymax = vi2->y; 
else 
if ( vi3->y > ymax ) 
ymax = vi3->y; 
LineSegment edge(*vi2, *vi3); 
edgesl.push_back (edge) ; 
vi2 = vi3; 
}// for i 
++pi0; 
}// while 
return ( edgesl); 


} 


// Intialises the Global Edge Table with the edges corresponding to the scan line y 
void SLAM::InitGET(GETable &get , vector< LineSegment> &edgesVector, float y) { 
GET tempGET; 
LineSegment e; 
float yinf, ymin, x; 
yinf = 1024; 
for ( unsigned int i =0; i < edgesVector.size(); i ++ ){ 
e = edgesVector[i]; 
ymin = min (e.vertex0.y, e.vertexl.y); 
if ( ymin == e.vertex0.y ) 
xX = e.vertex0.x; 
else 
xX = e.vertexl.x; 
if (ymin < yinf) 
yinf = ymin; 
if (yinf < 0) 
yinf = 0; 
y = yinf; 
tempGET.11 = e; 
tempGET.slope = (e.vertex0.y - e.vertexl.y ) / ( e.vertex0.x - e.vertexl.x ); 
tempGET.xOfymin = x; 
tempGET.ymax = max ( e.vertex0.y, e.vertexl.y); 
tempGET.ymin = ymin; 
if ( tempGET.slope != 0 
get.push_back(tempGET ); 
} 
sort ( get.begin(), get.end()); 
} 


// Updates the AEt entry with new edges depending on the scan line y 
void SLAM: :UpdateAET(AET &aet, float y, GETable é&get ){ 
AETEntry tempAET; 
aet_it ae = aet.begin(); 
int size = aet.size(); 
if ( !aet.empty()) 
for (int i = 0; i < size;itt+ ) 
{ 
if( ae->ymax < ytl 
aet.erase (ae); 
else 
aett; 
} 
get_it ge= get.begin() 
for ( unsigned int j= 


{ 


0; 3 < get.size(); 43 ++ ) 


if ( ge->ymin == y ) 

{ 
tempAET.edge = ge->11; 
tempAET.xval = ge->xOfymin; 
tempAET. ymax = ge->ymax; 
tempAET.m_inverse = 1/ge->slope; 
aet.push_back ( tempAET ); 

} 


gett; 


fis 


} 


// updates the x-val in the AETEntry with scan conversion 
void SLAM: :UpdateAETEntry( AETEntry &g ){ 


if (( g.xval > max( g.edge.vertex0.x,g.edge.vertexl.x)) ||( g.xval < min( 
g.edge.vertex0.x,g.edge.vertexl.x) )) 
return; 
g.xval = (g.xval + g.m_inverse) ; 


} 


// updates with Global Edge Table 
// removes the edges with ymin less than the scan line 
void SLAM: :UpdateGET( GETable &get, float &y) { 
get_it ge = get.begin(); 
int size = get.size(); 
for (int i = 0; i < size;itt+ ) 
{ 
if( ge->ymin <= y ) 
get.erase (ge); 
else 
gett; 
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